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Abstract 

Many densely packed suspensions and colloids exhibit a behavior known as Discontinuous Shear 
Thickening in which the shear stress jumps dramatically and reversibly as the shear rate is increased. 
We performed rheometry and video microscopy measurements on a variety of suspensions to determine 
the mechanism for this behavior. We distinguish Discontinuous Shear Thickening from inertial effects 
by showing that the latter are characterized by a Reynolds number but are only found for lower packing 
fractions and higher shear rates than the former. Shear profiles and normal stress measurements indicate 
that, in the shear thickening regime, stresses are transmitted through frictional rather than viscous 
interactions, and come to the surprising conclusion that for concentrated suspensions such as cornstarch 
in water which exhibit the phenomenon of Discontinuous Shear Thickening, the local constitutive relation 
between stress and shear rate is not necessarily shear thickening. If the suspended particles are heavy 
enough to settle we find the onset stress of shear thickening r min corresponds to a hydrostatic pressure 
from the weight of the particle packing where neighboring particles begin to shear relative to each other. 
Above Tmi n , dilation is seen to cause particles to penetrate the liquid-air interface of the sheared sample. 
The upper stress boundary T max of the shear thickening regime is shown to roughly match the ratio of 
surface tension divided by a radius of curvature on the order of the particle size. These results suggest a 
new model in which the increased dissipation in the shear thickening regime comes from frictional stresses 
that emerge as dilation is frustrated by a confining stress from surface tension at the liquid-air interface. 
We generalize this shear thickening mechanism to other sources of a confining stress by showing that, 
when instead the suspensions are confined by solid walls and have no liquid-air interface, T max is set 
by the stiffness of the most compliant boundary which frustrates dilation. All of this rheology can be 
described by a non-local constitutive relation in which the local relation between stress and shear rate 
is shear thinning, but where the stress increase comes from a normal stress term which depends on the 
global dilation. 

1 Introduction 

Shear thickening is a category of non-Newtonian fluid behavior in which the viscosity <r\ = r/ j increases as a function 
of shear rate 7 or shear stress r over some parameter range. A particularly dramatic manifest ation characterized by a 
sharp jump in s tress with increasing shear rate, often called Discont inuous Shear Thickening (|Metzner and Whitlockl 
1958; Hoffman. [19721 ; iBarnesI , 1 19891 ; llVIar anzano and Wagnerl . l2001al ). occurs in many densely packed suspensions and 
colloids such as cornstarch in water. These suspensions feel like a thin liquid at low stresses, but become very thick 
and can even crack like a solid at higher stresses, and become thin again when the stress is removed. Such materials 
are of practical interest for their properties as dampeners and shock absorbers (|Lee et all 120031 ; IShenov et all 12005 ; 



Jollv and Bender, 2006). While some milder types of shear thickening can be explained as viscous (|Bradv and Bossi: 



1985; Wag ner and Brady . 2009; Ch eng et alll201lf ) or inertial (|Bagn old. 1954) effects, prior approaches have not been 
very successful at describing the dramatic effects of Discontinuous Shear Thickening. Our goal with this paper is to 
use observations to develop and test a model for Discontinuous Shear Thickening. This problem can be broken down 
into several questions: under what conditions will a suspension exhibit Discontinuous Shear Thickening?, what are 
the scaling laws that determine the parameter range of shear thickening?, and what is the form of the constitutive 
law? We will use a wide variety of rheometry and video microscopy measurements to answer the above questions. 
Our approach differs from previous work in that we consider a granular point of view in addition to the traditional 
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hydrodynamic approach by investigating the source of observed stresses during Discontinuous Shear Thickening, 
rather than focusing only on the viscous contribution to the local constitutive relation. 

The remainder of this paper is organized as follows. In Sec. [5] we review the literature on shear thickening to 
characterize Discontinuous Shear Thickening and summarize existing models. In Sec. [3] we describe the rheometry 
techniques and suspensions used in our experiments. In Sec. 0] we show viscosity curves for shear thickening sus- 
pensions at different packing fractions and liquid viscosities. With this we can test inertial hydrodynamic scalings 
which we find to apply in a different parameter regime than Discontinuous Shear Thickening. In Sec. [5] we show 
measurements of the onset stress T m in of shear thickening for different particle sizes and liquid densities and show 
that for settling suspensions, the onset of shear thickening occurs when the shear stress becomes large enough to 
initiate relative shear between particles, balancing against a gravitational pressure. In Sec. [6] we show shear profile 
measurements with levels of density matching and demonstrate that the large stress jump that characterizes Discon- 
tinuous Shear Thickening is not directly affected by the local shear rate or gravity-induced inhomogeneities. In Sec.[7| 
we show shear and normal stress measurements with different boundary conditions. We show that the shear stress 
is coupled to the normal stress as in a frictional constitutive relation, rather than viscous, and that the rheology is 
controlled by the normal stress boundary condition. Combining this with the shear profile measurements, we come to 
the surprising conclusion that for concentrated suspensions such as cornstarch in water which exhibit the phenomenon 
of Discontinuous Shear Thickening, the local constitutive relation between stress and shear rate is not actually shear 
thickening. In Sec. [8] we show images and movies of a visible change of the suspension-air boundary as particles 
penetrate the liquid-air interface in response to dilation of the granular packing under shear, which coincides with the 
shear thickening regime. In Sec. H] we propose that this change in boundary condition due to dilation results in a con- 
fining stress from surface tension at the liquid-air interface, providing the normal stress boundary condition required 
for stress jump that characterizes Discontinuous Shear Thickening. We then show that comparisons of the measured 
stress and dilation are consistent with this model, and that the confining stress scale from surface tension agrees 
with measurements of the maximum stress in the shear thickening regime T ma x for a wide variety of suspensions. 
In Sec. [10] we generalize this result to other sources of confining stress with measurements from a rheometer with 
solid walls instead of a liquid-air interface at the boundary and show that T max is generally set by the stiffness of the 
boundary. Finally, in Sec. [Tljwe discuss a generalization of the stress scales that bound the shear thickening regime 
to parameter regimes where other forces are relevant, the general constitutive relation, connections to the physics of 
other types of materials, and summarize the conditions for the occurrence of Discontinuous Shear Thickening. 



2 Background 



Because of the vast amounts of literature referring to different phenomena and mechanisms which all fall under the 
category of shear thickening, we first carefully define what we mean by Discontinuous Shear Thickening. We will 
suggest based on a literature review that hydrodynamic models which have been successful for describing other types 
of shear thickening have been insufficient for describing Discontinuous Shear Thickening. Rather, there is significant 
evidence suggesting that dilation and normal forces play an important role which leads us to look for sources of stress 

similar to those in granular systems. 

Discontinuous Shear Thickening dMetzner and Whitlocl^ll958l;lHoffmanLll97lll974l.[l^;lBarnesl.[l989l;lBoersma et all 



1990; iLaunl Il994l ; iFrith et all . Il996 
h 11 l l„ , , , 'i V 



Bender and 



Wagnerl 

I onno 



1996 
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2000; Maranzano and Wagner, 
Shcnov and Wagner, 2005: 
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2001allbl. 120021: iBertrand et al.l.l2002l ; Lootens et al.l. 12003. 120051; lEgres and Wagnerl, 120051 
Egres et all . 120061 : iLee and Wagnerl l200d : iFall et all 120081 ; iBrown and Jaegerl . 120091 : iBrown et al, 
characterized by a set of several properties, which are observed to be similar in both dense suspensions and colloids. 
Thus we find it convenient to define the phenomenon of Discontinuous Shear Thickening based on these properties 
and to distinguish it from other types of shear thickening: 



1. Stress scales : The boundaries of the shear thickening regime are simply de scribed in terms of stres s scales 



(rather than shear rate) which are mostly independ ent of packing fractio n dMaranzano and Wagner]. [20 01a 
She nov and Wagnerl . 120051 : IBrown and Jaegerl . [2009) and liquid viscosity (|Boersma et all Il990l : IFrith et al. 



19961 ). Consequently, the onset shear rat e varies with pack i ng fra ction and liquid viscosity since suspension 



viscosities increase with both parameters (|Bradv and Bosslsl . [l985). Shear thinning or Newtonian behavior is 
found at stresses below the onset of shear thickening at a stress Tmin- There is an upper bound of the shear 
thickening regime at a stress T max , and shear thinning is usually found at higher stresses. 

Diverging slope: The term 'Discontinuous' refers to the apparent discontinuous jump in the stress r(-y) of 
orders of magnitude in the shear thickening regime. This jump is only observed at very high particle packing 
fractions <j) (around 0.5 for nearly spherical particles). The slope of r(7) is only very steep over a small range 
in packing fraction (a few percent), lessening significantly at lower packing fractions (|Maranzano and Wagnerl . 
l2001al : lEgres and Wagnerl 120051 : [b rown and Jaegerl . 120091 ) . The packing fraction dependence of the slope of 
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r(-y) can be characterized as a power law diverging at a packing fraction (j> c (|Brown and Jaegerl . [2009), For 
4> > 4> c the system is jammed meaning it will not flow for applied stresses below a non-zero yield stress of scale 

3. Reversibility: The Discontinuous Shear Thickening described above is reversible, meaning viscosity curves 
are similar whether they are measured with increasing or decreasing stress histories. Some examples of dra- 
matic shear thickening have been foun d to be irreversible because of chemi cal-attraction-ind uced aggregation 
(|Osuii et al.l . 120081 : lLarsen et al.l . [2010) or occur only in transient behavior jFall et all |2010| ). and will not be 
considered here as they may be different phenomena. 

Because the above features of Discontinuous Shear Thickening are similar in both the suspension and colloid 
regimes, we refer to data in both regimes even though the dominant forces may be different. 

To facilitate an understanding of Discontinuous Shear Thickening, it may be useful up front to appreciate a dis- 
connect between local and global viewpoints of rheology. Standard rheology experiments measure a global mechanical 
response based on the drag force required to move two solid surfaces at some speed relative to each other with a 
fluid in between. This is in contrast to the hydrodynamic theory of rheology which is based on continuum equations 
consisting of local constitutive relations between shear stress and shear rate. One of the surprising results presented 
in this paper is that the global and local viewpoints of rheology lead to drastically different stress vs. shear rate 
relations in the Discontinuously Shear Thickening suspensions being investigated. Historically, Discontinuous Shear 
Thickening has been reported in experiments based on measurements of the global mechanical response. We will show 
in Sections [5] that the constitutive relation based only on the local shear rate can actually be shear thinning instead 
of shear thickening for these suspensions. To resolve this apparent contradiction, in Sec. [7] we will show that most 
of the shear stress is due to frictional interactions which depend on the normal stress and can be controlled by the 
boundary conditions rather than the local shear rate. Because the two standard viewpoints disagree on whether these 
dense suspensions are shear thickening or shear thinning, there will be some difficulty with terminology. We chose to 
continue using the terminology of Discontinuous Shear Thickening to keep the connection to the previous literature, 
where the term is used extensively for similar experimental results. We have chosen to capitalize 'Discontinuous Shear 
Thickening' to refer to it as a name for a phenomenon rather than a local description, since from the local viewpoint 
it would not qualify as shear thickening. Towards our goal of understanding Discontinuous Shear Thickening, our 
approach is to investigate all sources of the measured stresses, rather than just the hydrodynamic contribution to 
local constitutive relations. 

In contrast to the Discontinuous Shear Thickening described above, there are other types of shear thickening 
with different characteristics that have been modeled in hydrodynamic terms. In the hydrocluster model, shear 
thickeni ng occurs when an increasing shear rate leads to particle clustering and consequently increasing viscous 
stresses (|Bradv and Bosslsl . Fl985; Wag ner and Bradvl . 120091 ). The shear thickening found in this model is relatively 
very weak in the sense that the viscosity increases by only a few percent per decade of shear rate |Bradv and Bossia . 
1985; Melrose and Balll . l2001al ). Inertial effects can also result in shear thickening, described in terms of either a 
Reynolds n umber or Bagn old number, in which inertial stresses can be characterized by T(j) oc j 2 in the limit of high 
shear rate (iBagnoldl . Il954l ). In each of these cases, the onset of shear thickening is characterized by a shear rate, in 
contrast to an onset stress for Discontinuous Shear Thickening. Because of the less steep r( , y), such shear thickening 
is often called Continuous Shear Thickening. Continuous Shear Thickening can be observed even at very low packing 
fractions and has a much weaker packing fraction dependence than Discontinuous Shear Thickening. In Sec. U we 
will show an example of inertial shear thickening and characterize the parameter regime where it occurs in so we can 
clearly separate Discontinuous Shear Thickening from inertial stresses. 

Because of the qualitative differences between Discontinuous and Continuous Shear Thickening, we suspect they 
are different phenomena with different mechanisms, and focus here only on Discontinuous Shear ThickeningQ For this 
reason, we will not assume that results that apply to Continuous Shear Thickening also apply to Discontinuous Shear 
Thickening. This is a different approach from much of the liter ature which has attempted to apply the hydrocluster 
model to both Continuous and Discontinuous Shear Thickening (|Bradv and Bossisl. 19851; O'Brien and Mackav , |2000| ; 



M aranzano and Wagned . [2001al ; lGopalikrishnan and Zukoskl 12004 lOsuii et alll200Sl ; IWagner and BradvH2009l ). The 



idea with the earlier approach has been to sta rt with the hydrodynamic models of Continuous Shear Thickening which 
are well-understood at low packing fractions (|Bergenholtz et all 120021 ) . and extend them to higher packing fractions. 
The expectation is that at higher packing fractions, hydrodynamically-induced clusters of particles form in which 
nearby particles act transientl y as a solid cluste r when they get so clos e to each other that the lubrication drag 
force between them blows up l|Bradv and BossiH . Il985l ; iFarr et all . 1 1997h . While this mechanism seems plausible, 



1 In some of the literature, suspensions and colloids which exhibit discontinuous shear thickening at high packing fractions 
are said to exhibit continuous shear thickening at lower packing fractions because the slope of stress vs. shear rate is lower. 
However, there is usually no qualitative change in behavior or identifiable transition in scaling when the packing fraction is 
varied. Thus we do not follow this convention because it does not suggest different phenomena and refer to such systems as 
exhibiting Discontinuous Shear Thickening. 
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so far the calculations have failed to reproduce steep viscosity curves comparable to experimental measurements of 
Discontinuous Shear Thickening. 

The major success of the hydrocluster model for Discontinuous Shear Thickening is the calculation of the stress at 
the onset of shear thickening r m i„. In early models for Brownian- motion dominated colloids the onset was described 
by a critical Peclet number Pe — Gnr/'ja 3 /kT for a particle size a and thermal energy k T. Shear thickening was 
expected to occur for Pe 3> 1 as the shear stress overcomes thermal diffusion of the particles jBradv an d Bossis. ll985l ; 
iMelrose and Ball l2001al ). However, observations found the onset of shear thickening to be determined by the same 
stress at different packing fractions rather than the same shear rate. This model, when using the suspension viscosity 
to convert to a stress scale T m in = rfj = fcT/67ra 3 , has b een successful at calculating the onset of both Continuous and 
Discon tinuous Shear Thickening in Brownian colloids (|Gopalikrishnan and Zukoskil . 120041 : iMaranzano and Wagnerl . 
2001b]). For colloids where electrosta t ic rep ulsions from a zeta potential £ are dominant the above model had to 
be modified (|Maranzano and Wagnerl . l2001aT ). In that case, the particular scaling found corresponds to a stress 
T m in ~ e^ 2 /" 2 for permittivity e which characterizes the electrostatic interactions between neighboring particlesQ 
In each regime, the modifications to the hydrodynamic model required to fit it to the data resulted in completely 
eliminating any dependence on hydrodynamic parameters such as viscosity or shear rate. The resulting onset stress 
scale is not restricted to hydrodynamic mechanisms, as any type of force transferred through a continuum system 
can be expressed in terms of a stress. Not surprisingly, with several relevant forces in colloids and suspensions, 
each of which could be dominant in different cases, a variety of different scalings for the ons et stress have been 
found. Depending on the parameter r ange, this domi nant force could be Brownian motion (|Bergenholtz et all 
|2002|; Gopalikrishn an and Zuk oski. 2004|), zeta potential (Maranzano and Wagnerl . [2001a! ). induced dipole attractions 



( Brown et al.1 . 2010al ) , or steric repulsion (|Hoffmanl . Il998 ) . Notably, in each case hydrodynamic terms such as shear 



rate and viscosity were found to be absent from the modified scalings required to match the experiments. This 
suggests inertia or hydrodynamics-based models are not necessary to determine the onset of Discontinuous Shear 
Thickening as initially envisioned by the hydrocluster model. 

Another observation that the hydrodynamic models require modification to describe is the upper stress boundary 
of the shear thickening regime T max , since for inertial effects the viscosity increases monotonic ally with shear rate . 
It has been suggested that this could be fixed by accounting for the finite stiffness of particles |Kalman et all [2009). 
but again this introduces a stress scale that is not necessarily a hydrodynamic in origin. We will address the issue of 
the upper stress bound r ma x in Sec. [5] 

We could simpl y approach the problem in terms of stress scales as a modification of hydrocluster models as many 
others have done (jMaranzano and Wagner] , l2001al ; lO'Brien and Mackavl . |2000| ; iGopalikrishnan and Zukoski |2004| ; 
lOsuii et all |200S| ). Instead, here we focus on the scaling laws for the stress scales r m i n and T max to gain insight 
into the relevant physical mechanisms and come up with a description that encompasses all of these scalings without 
the need to refer to a hydrodynamic model. A major advantage of this approach is that different mechanisms such 
as interparticle forces, gravity, and surface tension can be simply expressed in terms of stress scales which can be 
compared to the measured stresses without the need to reference a base model. Another good reason to try to 
understand Discontinuous Shear Thickening without reference to hydrodynamic models is that some features of the 
rheology suggest a granular rather than hydrodynamic mechanism. While a granular point of view is not necessarily 
the only way to understand the phenomenon, it will allow us to easily interpret many features. 

Granular materials can have properties of solids, liquids, or gases under different conditions l| Jaeger et all fl996h . 
For example, randomly packed particles at high enough packing fractions cannot shear or compress without deforming 
particles because geometric constraints force them to be in contact, so they have a yield stress like a solid. At lower 
packing fraction the particles are able to move around each other fre ely in a liquid-like sta t e. The transit i on be tween 
these two regimes is sharp and is known as the jamming transition (|Liu and Nage l. 1998; O'He rrTet all [2003). For 
shear thickening suspensions, the div ergence of the slope of r j-y) at a critical packing fraction (f> c was found to 
correspond to the jamming transition (|Brown and Jaeger] . 12009] ). 

Forces tend to be transmitted through jammed granular packings along concentrated paths called force chains 
dCates et all 19981) such that the distribution of for ces is characterized by an exponential tail (|Mueth et all 1 19981 ; 



iMaimudar and Behringerl 120051 ; ICorwin et all 120051 ). Simulations of shear thickening colloids have similarl y foun d 
contact networks between particles with an exponential distribution of forces under shear (Melrose and Ball, 2004b). 

A feature of granular shear flows with special relevance to shear thickening is dilation (|Revnoldsl ll885: On oda and Linigerl . 

1990). When a granular packing is sheared, the particles have to move around each other so the packing dilates, 
taking up more volu me than it does at rest. It has long been known that dilation occurs along with Discontinuos 
Shear Thickening (see lMetzner and Whitlo"ckl (1 19581 ). an d references th erein). In fact, in some of the literature 'dila- 
tancy' has been used as a synonym for shear thickening l|Barnesl . [l989l) . Especially important in understanding this 



2 While the forces were calculated using a hydrodynamic model with particles interacting by a lubrication force, the hydro- 
dynamic terms cancel up to a dimensionless coefficient of order 1. Since the model was an order-of- magnitude calculation, it 
would have resulted in just as good a fit with the data if the interactions forces were calculated for a different particle separation 
distance or with a different mediating mechanism. 
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relationship was the paper bv lMetzner and Whitlockl . They showed that for suspensions of 0.2-1 Ti02 particles, 
dilation initiated at stresses close to the onset of shear thickening for a range of packing fractions. However, they 
found suspensions of 28-100 /im glass spheres in sucrose solutions dilated but did not shear thicken, showing that 
dilation was not always equivalent to shear thickening. For 40+ years following this result, many of t he major papers 
on shear thickening dropped the focus on dilation as a mechanism in favor of viscous mechanisms ( Hoff man! . 1 19821 ; 
iBradv and Bossisl . 1 19851 : [Maranzano and Wagnerl . l2001ah . 

However, there is another possible interpretation of the data presented bv lMetzner and Whitlockl (|l958l ). They 
confirmed for several suspensions that the onset of shear thicken ing coincided with dilation, as had been seen in many 
previous results (see references in lMetzner and Whitlockl (j 1958)). Taken together with the observation of dilation in 
the absence of shear thickening, this inductively suggests that dilation is necessary but not sufficient for Discontinuous 
Shear Thickening. More recent results h ave shown th at shear thickening can be hidden by a yield stress or other 
shear thinning effect |Brown et al.l . l2010al ). Specifically. iMetzner and" Whitlock used glass beads ranging from 28-100 
fim in diameter in a Couette geometry, a nd such large, he avy particles will jam in a Couette cell because they settle 
under gravity, resulting in a yield stress ([Fall et al.. 2009). This can explain why shear thickening was not observed 
for the settling particles used by Me tzner and Whitlockl . I t has only been in the last 10 years that dilation has 
become p r omine nt again in the shear thickening literature (|0'Brien and Mackavl. l200d : lLootens et all . l2003t 120051 : 
iFall et all . 120081 ). However, a mechanism by which dilation leads to the dramatic increase in stress in Discontinuous 
Shear Thickening has yet to be explained. This is the subject of Sections [8l [9] and 1 101 

When dilation of granula r shea r flows is prevented by confi nement, shear is instead accompanied by normal forces 
against the walls (|Revnoldsl Il885l : lOnoda and Linigerl . Il990l ). Dilation of packings against boundaries can play a 
dominant role in the mechanics of granular system s where confining pressures from the boundary are transmitted 
through the material |Lambe and Whitman! . 1 19691 ). In Discontinuous Shear Thickening suspensions normal forces 
are usually found to be positive, meaning the sample is pushing against the rheometer plate a s expected for dilation 
dJomha and Reynoldsl 1 19931 : lLootens et all . 120031 . 120051 : IFall et all |200S| ; iBrown et all l2010bl ). It was proposed by 
iFalT et al.l (2008) that shear thickening cannot occur if the normal stress is taken away. We will test this in Sec. [7] 



3 Materials and methods 
3.1 Suspensions 

We studied a wide variety of suspensions with different particle sizes and shapes, liquid viscosities, and density 
differences among other properties to investigate the mechanism for Discontinuous Shear Thickening. We use this 
variety as a way to determine which features are common to all of the suspensions that exhibit Discontinuous Shear 
Thickening. 

As a prototypical shear thickener we used cornstarch obtained from Argo. Cornstarch particles have a mean 
diameter of 14 (jm and density of 1.59 g/mL based on buoyancy in CsCl solutions. They are very hydrophilic and 
hard - with a compression modulus on the order of 10 10 Pa - at room temperature. At higher temperatures the 
polymers that compose cornstarch particles can gel. To compare suspensions with different liquid viscosities, we 
suspended cornstarch in either a mixture of 61.5% water and 38.5% CsCl by weight with a viscosity of 1 mPa-s and 
density of 1.41 g/mL or a mixture of 73.5% glycerol 13.0% water, and 13.5% CsCl by weight with a viscosity of 80 
mPa-s and density of 1.34 g/mL. For such small particles, the settling time is several hours even without density 
matching. 

For a series of suspensions in which we varied particle size we used soda-lime glass spheres with a density of 2.46 
g/mL. We obtained particles with nominal diameter ranges of 3-10 /mi, 10-25 /im, and 15-40 /im from Corpuscular, 
45-63 /im, 75-104 fim (referred to as 100 /im), 177-250 (im, and 400-595 /im from MoSci (Class IV), and 1120-1350 
fim and 1900-2100 /im. For p olydisperse suspensions, measure ments of T m ; n could be collapsed onto the same curve 
as monodisperse suspensions (|Maranzano and Wagnerl . l2001al ). suggesting that polydispersity does not need to be 
accounted for, so we will compare particle distributions only by their mean diameter. The glass spheres were dispersed 
in various liquids, including water or mineral oil with a viscosity of 58 mPa-s and density of 0.87 g/mL. 

For visualization purposes we used two types of opaque particles. The first, referred to as Zr02, were spheres 
obtained from Glen Mills consisting of 69% Zr02 and 31% Si02- They have a nominal diameter range of 100-200 
/xm and a density of 3.8g/mL. These particles were dispersed in the same mineral oil used for the glass particles. 
For experiments with density matched suspensions, we used polyethylene spheres obtained from Cospheric. These 
particles have a nominal diameter range of 125-150 fim and density of 1.01 g/mL. They were dispersed in silicone 
oil AR 20 with a nominal density of 1.01 g/mL and viscosity of 20 mPa-s. When varying the temperature of the 
suspension, we found that the settling time was minimized at 19° C from which we estimated a density difference of 
order 10 -4 g/mL. 

We also include some summary data for particles fabricated into different shapes from polyethylene glycol (PEG) 
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suspended in liquid PEG-250 from iBrown et all (|201ll l. and 100 //m spheres made of polystyrene dimethyl ether 
suspended in PEG. 



3.2 Tools 

Measurements were performed with an Anton Paar Physica MCR 301 rheometer which measures the torque T required 
to shear a sample at a tool angular rotation rate u. Most measurements were done in a parallel plate setup where 
normal forces could be measured. This geometry is shown in Fig.QJ, and characterized by the plate radius R = 12.5 or 
25 mm and gap size d between the plates. A few measurements were done with a cylindrical cup-and-bob (Couette) 
geometry in which the environment is better controlled. The tool surfaces are smooth stainless steel. The viscosity, 
indicating the mechanical resistance to shear, is defined as r\ = r/7 in a steady state. For the parallel plate setup, 
for example, we represent the global shear stress by 



and shear rate by 



These standard coefficients in equations [1] and [2] correspond to the values at the outer radius of the plate in the case 
of a Newtonian shear profile. Since the dramatic feature of Discontinuous Shear Thickening is the increase in stresses 
under shear, the reported global shear stress and shear rate values are meant to characterize this response, in which 
the viscosity is a measure of mechanical energy dissipation, rather than a local constitutive relation between shear 
stress and shear rate. We will show in Sec. |6] that the local and global constitutive relations are quite different for 
Discontinuous Shear Thickening suspensions. The globally averaged stress is still appropriate to characterize forces 
in a way that is independent of system size so that measurements can be compared to the scale of forces associated 
with different physical mechanisms. Since most of the suspensions used are not density matched and non-Brownian, 
they can be very inhomogeneous. In this case, the given stress and shear rate relations may overestimate the average 
values by as much as 25% depending on the inhomogenieties. Since the formalism of local hydrodynamic constitutive 
relations and normal stress differences assumes homogeniety, we will not assume it applies, and instead base our 
interpretations on the gross mechanical response. 

In the parallel plate measurements, the upward force on the rheometer tool is measured and the mean normal 
stress rjv is obtained by dividing this normal force by the plate cross-sectional area. The standard deviation of the 
force measured during a static measurement over 10 s intervals with or without sample is 6 x 10~ 4 N, giving an 
uncertainty of 0.3 Pa (1.2 Pa) for a 50 mm (25 mm) diameter plate within a single measurement run. The standard 
deviation of the average force measured after calibration with no sample is 4 x 10 -3 N, giving an uncertainty of 2 
Pa (8 Pa) for a 50 mm (25 mm) diameter plate when comparing different runs. The resolution of the shear stress is 
much better, with an absolute uncertainty less than (0.001 Pa) 0.01 Pa for a 50 mm (25 mm) plate. 

The Anton Paar MCR 301 rheometer has special settings for normal force control measurements. Our reported 
measurements were done with the value of 'normal force hysteresis' set to 0.001 N. This value controls how much the 
normal force can deviate from the set value before the plate moves in response, although in practice the plate tends 
not to move until the normal force deviation exceeds about 0.01 N (20 Pa for the 25 mm plate). The 'normal force 
dynamics' value was set to 0% (default). This value controls the acceleration of the gap size in the feedback loop. 
Under this setting the gap responds slowly to variations in the normal force, and produces the most reproducible 
steady state gap sizes. 

The typical parallel plate setup is shown in Fig. in which the suspension is held in place between the parallel 
plates by surface tension at the liquid-air interface around the side. Because we saw that particles could penetrate the 
liquid-air interface (see Sec. [8} which we suspected could modify the stresses on the suspension from the interfacial 
tension, we desired a different boundary condition with a hard wall for some experiments. To accomplish this we 
machined an aluminum cylindrical cup with inner diameter 25.5 mm as shown in Fig. [Tp. The cup fits around the 
tool with a gap small enough to prevent 500 /jm diameter particles from slipping through but large enough to allow 
the tool to rotate without friction. This cup confined the particles to the volume beneath the plate, while the liquid 
could be filled to a higher level or omitted altogether so there was no liquid-air interface for particles to penetrate. 
The plates were covered with sandpaper sheets with a grit size of 100 fim to avoid slip with dry grains. In some cases 
we inserted a soft layer in between the top plate and the sandpaper to modify the compliance of the wall. 



3.3 Measurement procedure 

Packing fractions <j) were calculated as the volume of solid particles over the total volume of particles plus liquid 
mixed together. The packing density in terms of the inverse of the available free volume per particle may decrease 
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Figure 1: (a) A standard parallel plate rheometer setup. The suspension is confined between the plates by 
the surface tension of the suspending liquid, (b) A modified parallel plate setup with solid walls around the 
sample. In this setup the use of a suspending liquid is optional, and the wall stiffness can be modified by 
inserting layers of different stiffness between the suspension and plate. 




Figure 2: (a)Apparent viscosity curves for 100 /im glass spheres in mineral oil taken with different measure- 
ment durations to show hysteresis loops. Open symbols: measurement duration of 40 s per decade of the 
stress ramp. Solid symbols: 500 s per decade. Up-pointing triangles correspond to increasing stress ramps, 
while down-pointing triangles correspond to decreasing stress ramps, (b) Characterization of the hysteresis 
as the geometric mean of the viscosity ratio between the decreasing and increasing ramps of the hysteresis 
loop, plotted for different ramp durations per decade of stress. Dashed line: a ratio of 1 between increasing 
and decreasing ramps corresponding to no hysteresis. Dotted line: ramp rate used for later steady state 
measurements 



slightly during measurements as the grain packing dilates. Above the jamming transition, this packing density may 
also be less than the measured packing fraction if air bubbles become trapped in the interior. Humidity also has a 
large effect on the amount of water adsorbed onto dry grains open to the atmosphere, especially cornstarch which 
is so hygroscopic that 10-20% of the weight of the 'dry' powder is from water. Comparison of density measurement 
techniques suggest that cornstarch is porous or that it may absorb CsCl, thus we report mass fractions <f> m for 
cornstarch rather than volumetric packing fractions. 

Measurements were made with the rheometer's bottom plate temperature controlled at 20° C. The room humidity 
ranged from 22% to 38%, although during individual experiments the humidity was constant. To minimize evaporation 
or adsorption of water from the atmosphere to the suspension during measurements, we used a solvent trap when 
the suspending liquid was water which enclosed the sample and a small amount of air around it by an extra layer of 
liquid. The enclosed air equilibrated with the sample to prevent further changes to the suspension. 

The gap size d for parallel plate measurements was us ually about 1 mm, l arge enough to avoid finite-size effects 
on the viscosity for particles around 100 /im in diameter ijBrown et alll2010bl ). We measured bulk shear thickening 
with both rough and smooth plates and did not find any difference in the shear thickening due to the plate surface. 
To directly measure slip, we used video microscopy to observe the shear profile at the outer edge of the plate. The 
results of these measurements are shown in Sec. |6] We visually confirmed for all reported measurements that the 
suspensions do not spill. Spillage was often the limiting factor in the maximum stress or shear rate applied for our 
measurements. 
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Figure 3: (a) Viscosity vs. stress curves for suspensions of cornstarch in water. Mass fractions cf> m arc 
shown in the key; higher curves correspond to larger <j) m , and <p c = 0.56. The solid line corresponds to a 
constant shear rate, (b) A rescaling of the data as r/r/mj vs. Reynolds number Re = pijd 2 /r] m which should 
result in a data collapse for hydrodynamic flows. The line corresponds to a scaling r ~ j 3 / 2 in the range 
100 < Re < 3000. 



Suspensions were first pre-sheared immediately before measurements for at least 100 seconds at shear rates above 
the shear thickening regime where the steady state flow is fully mobilized, then viscosity curves were measured by 
ramping the control parameter (shear stress or rate) down and then up to obtain hysteresis loops. To ensure that 
we obtain steady state viscosity curves, the measurement ramp should be long enough that the size of the hysteresis 
loop is equal to that of the infinite duration limit. To check this, we show data in Fig. [2] for a sample of 100 pm 
glass spheres in mineral oil at (f> = 0.56, which is a stable sample over long time periods because the oil does not 
evaporate. Viscosity curves shown were taken first with a decreasing stress ramp followed by an increasing stress 
ramp for several different ramp durations. Since the control is a logarithmic ramp in stress over 4 decades, the 
measurement duration is specified in terms of duration per decade of stress. The hysteresis effect is characterized by 
the average distance between the upper and lower branches of the hysteresis loop in tj(t) on a log-log scale. This is 
calculated equivalently as the geometric mean of viscosity ratio (7j-(7")/7/+(r)]) 9 where rj- and rj+ are viscosities for 
decreasing and increasing stress ramps, respectively, and (...) s indicates a geometric mean, i.e. averaged on a log-log 
scale. This average was done over the stress range of 0.1 to 8 Pa in the shear thickening regime. The hysteresis 
initially decreases with increasing measurement duration, then levels off for long measurements indicating a steady 
state limit. The initial duration-dependent behavior is characteristic of a transient relaxation, and the crossover 
between the regimes indicates a characteristic timescale for the sample to reach steady state. The leveling off of the 
viscosity ratio at a value of 1 suggests that in this case there is not a true hysteresis effect. At packing fractions very 
close to the jamming transition we sometimes find some non-zero hysteresis loop even for very long measurements, 
such that different steady states can be reached dependent on the shear history. The steady state viscosity curves 
we report are generally in the long-duration regime where the viscosity ratio has leveled off, for this sample we use a 
control ramp rate of 500 s per decade of stress. For each steady state measurement, we ramped the control parameter 
down then up at least once, and in some cases up to five times, but we show only one set of curves for brevity if they 
were all identical within typical variations of 10-20% from run to run. 



4 Inertial scalings 

In this section we show some examples of viscosity curves for suspensions at different packing fractions. While there 
are already many examples of the packing fraction dependence, here we vary the liquid viscosity and the packing 
fraction from near the jamming transition all the way down to zero to compare to hydrodynamic scalings that apply 
to suspensions at low concentrations. 

We first show measurements of cornstarch suspended in water. These measurements were made with the Couette 
geometry. Viscosity is plotted vs. shear stress for several packing fractions in Fig. [3^,. Apparent shear thickening 
is seen as regions with a positive slope of the viscosity curve for all packing fractions, including <j> m — which is 
pure water. For pure liquids this behavior can be quantified in terms of a dimensionless Reynolds number, which 
represents a ratio of inertial to viscous stresses. This Reynolds number is usually of the form pid 2 -y/r]i for pure 
liquids where pi and r\i are the density and dynamic viscosity of the liquid, respectively. Thus, for the pure liquid 
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Figure 4: Viscosity curves for cornstarch in a glycerol-water mixture in which the Reynolds number remains 
small. Mass fractions <p m are shown in the key; higher curves correspond to larger <f> m , and <j) c — 0.58. The 
solid line of slope 1 corresponds to a constant shear rate and the steepest possible steady state viscosity 
curve. The vertical dashed lines define the stress scales T m i„ and r max that bound the shear thickening 
regime. 



the transition from a viscous-dominated regime with a nearly constant viscosity to an inertia-dominated regime with 
apparent shear thickening occurs at a fixed shear rate. However, we find that the sharp transition does not occur at 
the same shear rate for different packing fractions; specifically for low packing fractions the onset shear rate increases 
with packing fraction. T he contribution of viscou s stresses to the viscosity of dense suspensions is much higher than 
that of the pure liquid (|Bradv and Bosslsl . 0-985), so the viscous term in the denominator of the Reynolds number 
should be modified for suspensions. For the contribution of viscous stresses to suspension viscosity we do not take 
the zero shear rate limit of the viscosity since in this limit suspensio ns rheo l ogy can be dominated by non-viscous 
particle interactions which result in shear thinning at low shear rates (|Barnesl . Il98 9 ; Mara nzano and Wagner, 2001a; 
iBrown et all l2010al h Rather, we take as our best estimate the minimum suspension viscosity r\ m which occurs at 
the onset of shear thickening for each packing fraction. Our suspension Reynolds number is then Re — pi^fd 2 /r/ m . 

In Fig. [3b we plot the shear stress normalized by »7m7 (corresponding to the viscous contribution to the stress) 
vs. suspension Reynolds number Re. With this non-dimensionalization, data for any Newtonian fluid should collapse 
onto the same curve, which is nearly flat at low Re in the viscous regime and transitions to a more positive slope 
at higher Re due to inertial effects. We find that the lower mass fractions 6 m ~ 0.4 indeed collapse onto a single 
curve. The very mild increase in effective viscosity for 1 < Re < 100 is typical of hydrodyn amic flows in this 
range of Re where v i scous stresses are dominant but inertial effects start to become measurable l|Schlichtind , ll960l ; 
iKulkarni and Morrisl . 120081 ) . The app roximate scaling r ~ 7 3 ^ 2 for 100 J; Re ~ 3000 as inertia becomes stronger 
is also typical of hydrodynamic flows jDonnellv and Simonl . Il960lh The asymp totic scaling r ~ pi 2 d 2 expected in 



the fully inertial regime is usually not found until Re > 10 l|Schlichtind.ll960h . The data collapse and scaling for 
4> m ~ 0.4 suggests that in this regime the suspension behaves like a Newtonian fluid. In contrast, for 4> m > 0.4 the 
normalized viscosity curves deviate significantly from this collapse and scaling, increasing more steeply than inertial 
or viscous stresses are expected to be able to grow, and the steep shear thickening onsets at much lower values of 
Re. The onset of shear thickening occurs at lower shear rates for larger A, instead of the higher shear rates which 
would be required before inertial stresses dominate over viscous stressesjj Instead, the onset appears to be set by 
a constant stress scale for <j> m > 0.4 as seen in Fig. 3a. This suggests that inertial stresses are relatively small in 
this regime where shear thickening becomes more dramatic. This does not necessarily imply that viscous forces are 
dominant in the system, rather the lack of data collapse by the hydrodynamic non-dimensionalization for <f> m > 0.4 
in Fig. suggests other stresses must be involved. 

The scaling seen in Fig.[3]suggests two competing mechanisms for different types of shear thickening - inertial and 
Discontinuous - such that only the stronger effect is observed in a single viscosity curve. To obtain a system where 
the Reynolds number remains low in the Discontinuous Shear Thickening stress range even at low packing fractions, 
a liquid of higher viscosity can be used. Accordingly, we suspended cornstarch in a glycerol-water mixture with a 



3 For (p m > 0.4, r\ m is likely an overestimate of the viscous contribution to viscosity as the non-Newtonian terms become 
larger near the jamm ing transition and r) m likely represents a cross-over between the shear thinning and shear thickening effects 
iBrown et a l.. 2010a) Regardless, since the viscous contribution to the viscosity increases while the onset shear rate decreases 
with packing fraction, the onset Reynolds number still becomes very low at high packing fractions. 
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viscosity 80 times that of water. The viscosity curves vs. stress for different packing fractions can be seen in Fig. [4] 
Comparing with the data for cornstarch in water in Fig. we can see that there is similar strong shear thickening at 
high packing fractions, but no apparent shear thickening at low packing fractions. At the same stress, the Reynolds 
number is lower by about a factor of the ratio of the viscosities squared (ss 600) for the suspension with glycerol, 
thus Re < 100 and inertial effects remain negligible in the stress range of Discontinuous Shear Thickening even in the 
limit of zero packing fraction. We can see the remaining shear thickening uncontaminated by inertial effects is now 
very weak at <j>m ~ 0.40, and is almost imperceptible at (f>m ~ 0.30. This is a typical example of Discontinuous Shear 
Thickening, in which the region with positive slope of T)(t) occurs in a stress range that is nearly independent of 
packing fraction. This slope increases with packing fraction, approaching r\ ~ r (solid line in Fig- [4j corresponding to 
a discontinuous stress/shear-rate relation. The bounds of the shear thickening regime are characterized on the lower 
end by r m j n defined as the onset of a positive slope of ??(r), and on the upper end by T ma x defined as the transition 
from positive to negative slope. These transitions are measured as the crossover between local power law fits on either 
side. Because of fluctuations typically on the scale of 10-20% in the viscosity we do not count any features smaller 
than that threshold as distinct transitions. 

There are several other dimensionless numbers that have been used to describe inertial effects in particulate flows. 
In particular, o ften the system size d is replaced with the particle size a, for example leading to a Bagnold number 
(Bagnold, 1954) or particle Reynolds number. The distinction between the two types of scalings can be made with 
a pure liquid at zero packing fraction which has the system size scale of d but no particle length scale. Since the 
cornstarch-in-water data for <f> m < 0.4, including <j) m — 0, collapse based on a Reynolds number scaling in terms of d, 
the system size should be taken as the relevant length scale as is typical in pure fluids, and a dimensionless inertial 
number based on particle size would not be able to collapse the data over as wide a range. 

We can qualitatively distinguish inertial flows from Discontinuous Shear Thickening because the packing fraction, 
shear rate, liquid viscosity, and gap size dependence of inertial flows all differ from that of Discontinu ous Shear 
Thick ening, and the steepest possible scaling for inertial flows is r ~ j 2 in the limit of large shear rates l|Bagnoldl . 
1954). Since the focus of this paper is on Discontinuous Shear Thickening, all of the following data will be in the high 
packing fraction regime (equivalent to <f> m > 0.4) and for Re < 100 to avoid mixing inertial effects with Discontinuous 
Shear Thickening. 

5 Gravity and the onset stress 

Suspensions and colloids occupy a region of phase space where many physical forces may be relevant; these include 
Brownian motion, gravity, surface tension, and electrostatics. One consequence of this is that there are different 
scaling laws for the onset of shear thickening in d ifferent parameter regimes where one of these forces is dominant 
llHoffmanl . 1 19981 ; iMaranzano and Wagneil l2001al lbl; iBergenholtz et all 12002] ; IShenov and Wagneit 120051 ; iBrown et all 
1201081 ). Here we address the case of a gravity-dominated regime for large particles. This regime is not yet as well- 
characterized as the other regimes, but an understanding of the effects of gravity on the onset stress will make it 
possible to further generalize the conditions for the onset of shear thickening with a mechanism that accommodates 
all of these scalings. 

We measured steady state viscosity curves for a series of suspensions of glass spheres (p = 2.46 g/mL) with 
diameters ranging from 6 to 2000 /im in diameter in either mineral oil (p — 0.88 g/mL) or water (p = 1.00 g/mL) 
at packing fractions ranging from 0.50 to 0.58 These particles are large enough to settle over time since the glass is 
much denser than the liquids. In each case, the major features were qualitatively similar to Fig. [4] with increasingly 
steep slopes at higher packing fractions and shear thickening occurring in a relatively fixed stress range. We obtained 
mean values of r m i„, corresponding to the onset of a positive slope of t](t), for each suspension averaged over a range 
of packing fractions close to but below <f> c . The mean values of T m i n for each particle size are plotted in Fig. [5] 

For the largest particles with a > 500 pm, the suspensions would not remain confined between the rheometer 
plates with a vertical boundary because the particles are so heavy that they can no longer be confined by surface 
tension. Since this confinement is set by the interplay between gravitational and surface tension forces, it is no 
surprise that loss of confinement occurs for particles on the order of the capillary length, or ~ 1 mm. We could still 
make measurements with some sample extended outside the area between the plates to obtain the scales of r m ; n and 
Tmax ■ For the largest glass beads with a diameter of 2000 pm, we found no shear thickening regime. The significance 
of this maximum particle size for shear thickening will be addressed in Sec 111.21 

It can be seen in Fig.0that there are two distinct scaling regimes for T m i n which meet at a minimum near a particle 
size of 50 pm. For smaller particles with a < 50 pm that approach the colloidal regime, interparticle interactions 
from various sources including electrostatics and Brownian motion tend to become large relative to gravity and can 
affect the onset stress. One effect which is relevant here is a high particle- liquid interfacial tension which results 
in an effective attraction between particles which can form force chains that span the system and jam it. This in 
turn results in a yield stress and shear thinning even at low packing fractions which then can hide shear thickening 
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Figure 5: The stress at the onset of shear thickening T m i n for glass spheres of different diameters a in mineral 
oil (solid circles, Ap — 1.58 g/mL) or water (open circles, Ap = 1.46 g/mL). Open diamond: polyethylene in 
silicone oil (Ap rs 0.0001 g/mL). The solid line is the shear stress required to lift particles off the top layer 
of the packing against friction and gravity p e //Apga/15.3. Dashed line: representative curve for data where 
gravity is not the dominant interparticle interaction. Dotted line: bound above which larger particles did 
not exhibit any shear thickening regime. 

(jBrown et all l2010ah . Specifically for this measurement series, with glass beads 50 /im and larger, the particles will 
disperse well and shear thicken in either oil or water. However, the 6 pm glass particles are effectively hydrophilic. 
Consequently in oil they have a significant yield stress and shear thickening was not observed at all. 

To understand the scaling of T m in for the larger particles with a > 50 pm, we now analyze the effects of gravity 
on non-Brownian suspensions. In the limit of zero shear rate, gravity results in particles settling and resting on the 
bottom plate. The measured stress would come only from shear of the thin fluid layer on top of the settled particles. 
The drag force from the shear in the liquid layer above the particles can start to move the upper layer of particles 
if it exceeds the static frictional force between particles under gravity. In a parallel plate geometry, the horizontal 
cross-section has a uniform area so, to balance forces, the shear stress r must be on average independent of height. 
As an estimate for t he drag force on a part icle in the top layer, we use the drag force on a sphere sitting on a flat 
surface is 2.55-kto 2 (|Goldman et all 1 19671 ) . The frictional force on one of these particles is irp e ffApga 3 /6 for an 
effective static coefficient of friction p e ff, density difference Ap between the particles and liquid, and acceleration of 
gravity g. Since the particles are sitting on a pile of beads instead of a flat surface, the effective friction is enhanced 
by geometry because of the need for the spheres to rise over the particles in the layer below. To measure p e ff, we 
glued 1 mm glass beads in a monolayer to a glass slide. We then performed an inclined plane test with this system 
immersed in water, setting individual glass beads on top of the bead-covered slide and slowly tilting the slide until 
the loose beads started falling down. From this we obtained p e ff = 0.8L0.1. Balancing the drag and frictional forces 
gives the stress at the onset of shear between grains to be p e ffApga/15.3. This prediction is plotted in Fig. It is 
seen to match well the measured onset stress r m in for particles between 50 and 1125 pm in diameter. This confirms 
that the onset of shear thickening in the gravity-dominated regime is set by the stress required to initiate shear of 
the particles against gravity and friction. 

Since the onset scaling for large particles is set by gravity, this suggests T m m can be lowered by density matching. 
We tested this by measuring steady state viscosity curves for 100 pm polyethylene particles in silicone oil with a 
density difference of about 10 -4 g/mL. The mean value of r m in is shown in Fig. [S]by the open diamond. While Ap 
was reduced by a factor of 10 4 compared to the non-density matched case, the onset stress was only reduced by an 
order of magnitude. In this case the density-matched data fall onto a similar scaling as was found for the glass for 
a < 25 pm. In many cases for particles even as large as 100 pm, we found density matching can have no effect on 
the onset stress due to the significance of interparticle interactions. For cornstarch in water, density matching by 
adding CsCl to the water did not reduce the onset stress. For glass spheres in a heavy liquid p = 2.46 g/mL (Cargille 
labs inorganic salt series) we found no measurable decrease in the onset stress compared to mineral oil or water, and 
found shear thinning below the onset stress as opposed to the Newtonian scaling found for glass suspensions whose 



11 



onset is determined by the gravitational scaling ijBrown and Jaeeed . [20091 ; iBrown et all , l2010ah . This suggests that 
in each of these cases, the stress scale characterizing interparticle interactions which is dominant for smaller particles 
is very close to the onset stress if not the dominant factor. These results show that while density matching can lower 
the onset stress in the gravity-dominated regime, it cannot do so beyond the limits set by any other stress scales due 
to particle interactions. Thus we generally expect a larger effect of density matching for very large particles further 
into the gravity dominated regime. 

In this set of experiments with a parallel plate setup, gravity caused particles to settle and shear thickening 
required initiating shear so the onset stress T m in scaled like a hydrostatic pressure due to the weight of the top layer 
of particles only. In contrast, in a Couette cell with vertical walls a yie l d stre ss was found that scales with the same 
hydrostatic pressure in the suspension pushing on the walls jFall et al.l . l2009T ). We know that such a yield stress can 
hide shear thickening if it is larger than the stress from shear thickening mechanisms (|Brown et al.l . l2010aT ). This 
implies that the effect of the yield stress from gravity is shear-geometry dependent because of the directionality of 
gravity . The fact that thi s yield stress can move the onset stress and hide shear thickening if the yield stress is 
larger (|Brown et al.l . l2010al ) suggests it works in addition to stresses responsible for shear thickening, and there is no 
indication that the shear thickening mechanism itself is affected by gravity. We will revisit the issue of the significance 
of the scalings for the onset stress in Sec. 111.21 after a mechanism for shear thickening is identified. 



6 Shear profile 

In this section, we show shear profile measurements of both density matched and non-density-matched suspensions 
that exhibit Discontinuous Shear Thickening. The inhomogeneity due to gravity creates a shear gradient that allows 
us to separate out the contributions from viscous forces, gravity, and other forces to the constitutive relation between 
stress and shear rate in the shear thickening regime. 

To measure the shear profile we used a video camera with a bellows and magnifying lens to obtain a pixel size as 
small as 10 /im. The camera was placed next to the standard parallel plate rheometer setup and focused on the outer 
edge of the sample in the plane of the shear direction and shear gradient. While there is some distortion from looking 
through the curved liquid-air interface, we can track individual particle motions to measure the shear profile at the 
edge of the sample. Videos were taken for constant shear rate conditions after the steady state was reached. Steady 
state shear profiles obtained by using particle image velocimetry to obtain local particle velocities and averaging the 
velocities at each height. A small tilt of the camera caused a smoothing effect over about 4% in the depth. 

We first describe results for a settling suspension of 150 /Urxi ZrC>2 spheres in mineral oil at (f> — 0.53 with a gap 
d = 890 fim. This suspension is chosen for visualizations instead of glass because the particles are opaque. The raw 
particle motions under shear are shown in supplementary videos 1 and 2 for two different shear rates. Velocity profiles 
are shown in Fig. [6^ for a range of shear stresses. The corresponding global viscosity curve is shown in Fig. [6)3. Below 
T m in ~ 0.3 Pa, we found no measurable particle motion up to a resolution of 10 times the plate displacement. In 
this regime the particles remained settled due to gravity as expected based on the measurements of T m in in Sec. [5] 
Above T m i n , we found a narrow shear band near the moving top plate . This also agrees well with t he observations 
that the onset of shear thickening corresponds to the onset of dilation (|Metzner and Whitlockl . I1958T I , since shearing 
of the grains is what results in dilation. The width of the shear band increased as the stress was increased. Layering 
was clearly observed at higher shear rates, which results in the step-like shear profiles in Fig. [6^i. Effects of th is 
layering on the measured stress are only expected for smaller gaps, less than about 5 layers ( Brown et al.l . [2010bl ). 
We performed similar measurements with glass particles in mineral oil with a gap 12 particles wide as opposed to 6 
particles with the Zr02- Results were qualitatively similar to the ZrC>2 data, although layering was less prominent 
as expected for a finite-size effect, appearing only clearly in the top two layers. 

We next describe results for nearly density-matched 135 /im polyethylene spheres in silicone oil at <f> = 0.55 with 
a gap 850 /im wide. The lighting used to take the videos heated the silicone oil by several degrees, so there was a 
slight density difference of about Ap = —0.01 g/mL such that the particles were slightly buoyant, effectively reversing 
the direction of gravity. Because the smaller density difference moves the onset of shear thickening to very low shear 
rates, we did not obtain measurements of the shear profile below T m in ~ 0.01 Pa. Velocity profiles are shown in 
Fig. [T^l for a range of shear stresses. The corresponding global viscosity curve is shown in Fig. [7p. In the shear 
thickening regime, the velocity gradient in the bulk was relatively small, with a shear band at the bottom plate and 
a layered structure at the top plate. In this case, the shear band appeared at the bottom plate and the direction of 
curvature of the shear profile was reversed due to the inversion of gravity. The shear band widened at higher stresses, 
similar to the case for the settling ZrC>2. Interestingly, for all three suspensions the layering was most pronounced 
near the top moving plate despite the gravity inversion. 

These measurements also allow us to measure slip directly. When the particles were settled with the plate moving 
past, there was not even contact so the difference between plate and particle motions is technically not slip and so 
would not be expected to follow slip correction models which usually assume a linear bulk velocity profile. Settling 
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Figure 6: (a) Shear profiles at the plate edge for settling particles of ZrC>2 in mineral oil (Ap = 2.9 g/mL). 
The mean velocity v is normalized by the plate edge velocity v p , and the depth h is normalized by the gap 
d. Shear stress r for each profile shown in the key; higher curves correspond to larger r. Dashed line: upper 
bound of 10 -3 for a measurement at r < r m i„. Dotted black line: depth equal to 1 particle diameter. Solid 
lines: fits of Eq. [S]to the data for r < T max . Inset: same data on log-linear scale, (b) Local viscosity curves 
based on the local shear rate from the shear profile. Open circles: local viscosity in bulk region. Solid circles: 
local viscosity in the shear band near the top plate. Solid line: global viscosity curve. 




Figure 7: (a) Shear profiles at the plate edge for density matched polyethylene in silicone oil (Ap w —0.01 
g/mL). The mean velocity v is normalized by the plate edge velocity v p , and the depth h is normalized by 
the gap d. Shear stress r for each profile shown in the key; lower curves correspond to larger r. Dotted 
black line: depth equal to 1 particle diameter. Solid lines: fits of Eq. [5] to the data for r < T max with the 
substitution h/d — > 1 — h/d since the particles are lighter than the liquid, (b) Local viscosity curves based 
on the local shear rate from the shear profile. Open circles: local viscosity in bulk region. Solid circles: local 
viscosity in the shear band near the bottom plate. Solid line: global viscosity curve. 
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and slip effects can be distinguished by comparing to the density matched case shown in Fig. [7^, where the settling 
rate is much lower than the shear rate for all of the data shown. The difference between the speed of the top plate 
and neighboring particles in the more developed flow regimes is around 25%, roughly independent of shear stress. 
This does not change significantly at the boundary between shear thickening and shear thinning regimes, confirming 
that those rheological boundaries are not determined by slip. Since the goal of this paper is to understand the global 
response of Discontinuous Shear Thickening, we do not 'correct' for slip. Making a correction for slip would not 
significantly alter the shape of the viscosity curves nor move the regime boundaries in terms of stress because it only 
affects the shear rate, although it would slightly shift the magnitude of the viscosities reported. The lack of contact 
between the particles and plate with settling is not problematic in terms of the mechanical response because viscous 
interactions within the liquid transmit stress between them just as well as hard contacts. This is confirmed by our 
observation that switching from smooth to rough plates does not change the stress scales or whether Discontinuous 
Shear Thickening occurs. 



6.1 Local constitutive relations 



Here we use the shear profile to test constitutive relations in the shear thickening regime. Since the shear stress in a 
parallel plate geometry is independent of height, a local hydrodynamic constitutive relation r(ffi) dependent only on 
a local shear rate 7; would correspond to a linear velocity profile. To explain a non-linear steady state velocity profile, 
models have been introduced i n the past to account for fl uctuations in the local shear rate to an effective kinetic 
temperature dNott and Brady . 19941; Bocquet et all. 200lf) and the effect of th e local variation in packing fraction 
on the viscosity ( Nott and Brady . 19941 ; Bocquet et all l200ll ; iFall et ah . 2010T). In granu lar shear flows, the initial 



inhomogeneity is usually attributed to dilation near the moving plate ( Mueth et all . 2000h 



Here, since the curvature of the shear profiles in Figs. [6] and [7] changed with the direction of effective gravity 
for the particles, w e suggest that here gravity and friction are responsible for curvature in the shear profile, as in 
sedimenting flows l|Lenoble et al.l . 120051 ). Specifically, there can be frictional forces between particles due to the 
weight of the packing which increases with depth h into the sample relative to the top plate (for downward gravity) 
if the particles remain in contact via force chains. A non-linear shear profile could be the result of such an explicit 
height-dependence. The simplest form for a local stress relation that includes gravity is 



T = r;, y 7; + Tgh/d + r c (3) 

where r\ v is the viscous hydrodynamic contribution to the viscosity, the gravitational stress scale r g = /j, e ffApgd/15.3 
from Sec. [5] and r c represents any stresses that are independent of local shear rate and depth such as interparticle 
attractions. Rearranging gives the local shear rate 



7i = (t - t c - Tgh/d) /r] u . (4) 

This implies a critical depth h c /d = (t — t c )/t 9 at which the shear rate equals zero and beyond which there is no 
shearing of grains. This suggests the shear stress must exceed the sum of gravitational stress and interparticle stresses 
(included in r c ) on the first layer of par ticles (r a a/d) to shea r grains, in agreement with the condition for the onset 
of shear thickening shown in Fig. [51 and [Brown et all (|2010al l. respectively. The velocity profile can be obtained by 
integrating the local shear rate from Eq.fi] over depth. There are two solution regimes: 
if h c < d, then 




and if h c > d, then 




where the plate velocity v p = dj and the viscous stress scale is defined by t u = 77^7. The curvature of the velocity 
profile characterized by a quadratic term is set by the ratio of gravitational to viscous stresses t 3 /t v . The velocity 
profile becomes linear in the limit where this ratio goes to zero (h c > d) as expected. These profiles are concave up, 
and become more linear with increasing r in qualitative agreement with the data in Fig. [IJl The equations were 
written for the case where the effective gravity on the particles is downward. For the polyethylene data where the 
effective gravity is upward, we have to make the substitution h/d — > 1 — h/d which reverses the concavity. 

Because we are applying a continuum model to a system that is quantized due to layering, and fluctuations could 
smooth out mean shear profiles, this model will only be able to crudely approximate the slope and curvature of the 
shear profile. To test this model, we fit the function 
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Figure 8: Quadratic curvature 02 obtained from fit of Eqn.[7]to velocity profiles. The curvature is normalized 
by the model prediction from Eq. [6] with T m % n used as an estimate for the viscous stress r v . Data is fit 
for different normalized shear stresses r/r maa . for glass in mineral oil (solid triangles, Ap = 1.58 g/mL), 
polyethelyne in silicone oil (open circles, Ap = —0.01 g/mL), and ZrC>2 in mineral oil (solid circles, Ap = 2.9 
g/mL). The data collapse close to a value of 1 for r < T max suggests the curvature of the shear profile is 
due to the weight of the particles on deeper layers which is transferred via frictional contacts, and that the 
contribution of viscous stresses to the viscosity does not increase significantly in the shear thickening regime. 



v/Vp = ai ; h a 2 

a 



(h c - h) 



d 



(7) 



to the measured velocity profiles for each shear rate. Some of these fits are shown in Figs. [6k and [7k. For ZrC>2 and 
glass, it appears that h c < d for r < r ma x so we fix a\ = according to Eq. [5] The quadratic coefficient 02 can be 
compared to the prediction of Eqns. [5] and [6] with an estimate for the viscous stress t v . We showed in Sec. [5] that at 
the onset of shear thickening the viscous stress must be just enough to initiate shear so r v w r min . Since the shear 
rate increases slowly in the shear thickening regime for Discontinuous Shear Thickening, we will use this estimate for 
the entire shear thickening regime. The measured curvature 02 normalized by the predicted value T g /2T m in is plotted 
in Fig. [8] for each fit for ZrC>2, polyethylene, and glass. 

In the shear thickening regime (r < T max ), the data for all three density differences collapse onto a single curve 
with ai) {jgjlTmin) ~ 1- This value is in agreement with the model which confirms the curvature of the shear 
profile is set by a balance of gravity-induced friction and viscous interactions throughout the shear thickening regime. 
This balance implies the weight of the packing builds up in deeper layers, which requires force chains of solid particle 
contacts must extend from plate to plate, which is a common feature of granular systems. Since the shape of the shear 
profile depends on the specific force balance, the quadratic shear profile prediction is specific to settling suspensions 
within this model. Because the data collapse works reasonably well for systems of different sizes, it rules out the role 
of finite size effects in setting the velocity profile curvature. 

The data collapse of the shear profile curvature to a constant value in the shear thickening regime also suggests 
that the contribution of viscous stresses to the global viscosity is not increasing significantly in the shear thickening 
regime, and remains close to T m i n , in contrast to the expectations of hydrodynamic models for shear thickening. If 
the stress increase in the shear thickening regime was due to viscous forces proportional to shear rate, the reduced 
curvature would have to follow a slope of —1 in Fig. [8] in the range r < T ma x- Rather, the increase in stress in the 
shear thickening regime must be hidden in the uniform term r c due to other so-far-unidentified forces in Eqn. [3] 

For t > T m ax , the different curvature values do not collapse onto a single curve, suggesting either that the model 
fails in this regime or at least that r m i„ is no longer a good approximation of r v in this regime. Rather, the profiles 
appear to be closer to exponential (see inset of Fig. [Bk), similar to granular shear profiles of spherical particles 
(|Mueth et all [2000). This suggests that above T max the shear profile could be that of a fully granular system where 
there is no need for a contribution of viscous hydrodynamics. 

Curvatu re in the shear profile has been attributed to variations in the local packing fraction in some other 
experiments (|Fall et allfeoich . In a hydrodynamic model, small changes in packing fraction from dilation and viscous 
resuspension are significant because of the viscosity divergence with packing fraction as the viscous lubrication layer 
goes to zero at the jamming transition. However, our results on the curvature of the shear profile suggest stress 
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is transmitted mostly through frictional contacts rather than viscous interactions. In granular mechanics, fric tional 
contact forces can vary by less than about 30% with changes in packing fraction (|Lambe and Whitman! . 119691 . and 
this contribution would be small in comparison to the separation of local viscosities seen in Figs. [5p and[7p and the 
variations in curvature in Fig. [8] Another reason to suspect that viscous forces would not account for the measured 
stresses in the shear thickening regime comes from the magnitude of the me asured viscosity. Since viscosity values 
measured are up to 10 7 times the solvent viscosity l|Brown and Jaegerl . [2009). this would require subatomic gaps in 
a lubrication mode l. However, lubrication in mo lecular liquids breaks down at 2 molecular layers, below which the 
liquid is frictional (|Van Alsten fc Granickl |l988) . These issues suggest the stresses in the shear thickening regime 
must be explained by some non- viscous mechanism. 

Another mo del attributes curvature in the shear p r ofile t o a gradient in kinetic energy due to fluctuations in 
particle motion (|Nott and Bradvl. 1 19941 ; iBocquet et al.l . [2001). The contribution of this effect to the local stress 
gradient can be estimated as Vr ~ pvX/v ~ pj 2 d using the observation that the scale of the rms fluctuations in 



veloc ity are comparable to the mean flow velocity in granular flows with solid particle contacts (|Bocquet et al.l . 
120011 ).' This contribution to the stress gradient is at most on the order of 10 4 times the gravitational contribution 
to the stress gradient Vr ~ Apg even at the maximum stress in the shear thickening regime for the measurements 
shown in Figs.[6ji and[3i. Thus we expect the contribution of the kinetic energy to the shear profile to be negligible 
in the regime of our measurements. 

6.2 Localized viscosity curves 

Here we investigate the validity of local constitutive relations by plotting local viscosity curves. The local viscosity 
can be calculated as the ratio of the measured global stress and the local shear rate from the derivative of the velocity 
profile. To separate the bulk region from the shear band regions we use the mean slope of the velocity profile over 
different ranges of depth. For the profiles in Fig. [fji we use the range 0.25 < h/d < 0.8 for the bulk, and h/d < 0.3 
for the shear band at the top plate. For the profiles in Fig. [JJi we use the range 0.2 < h/d < 0.8 for the bulk, and 
h/d > 0.85 for the shear band at the bottom plate. According to the model of Eqns. [5]anc[6] the local bulk viscosities 
roughly approximate the linear term, which corresponds to the viscous contribution to the viscosity. These local 
viscosities are plotted in Fig. \5jp and [7p, respectively, along with the global viscosity curves. For each suspension, 
the shear band shows shear thickening similar to the global curve, while the bulk region appears to be everywhere 
shear thinning based on the local viscosity. In the non-density matched case below r m i„, the bulk was observed to 
be settled (i.e. locally jammed), corresponding to an infinite local viscosity. Thus the bulk region appears to have 
a yield stress due to gravity and consequently shear thinning but no shear thickening. The region that qualitatively 
determines the global rheology is not the bulk but rather the near-wall region where the shearing occurs. This is 
not surprising from a granular or solid mechanics point of view where the global behavior is often determined by the 
failure in the weakest region. 

Another test of the validity of local constitutive equations comes from comparisons of measurements in different 
measuring geometries. For a measurement in a cone and plate geometry the mean shear rate is independent of 
radius because the plate speed is proportional to the gap height at each point along a radius, while in a parallel 
plate setup the mean shear rate increases with radius because the plate speed is faster near the edge but the gap 
remains the same. Thus, assuming a local constitutive relation between shear stress and rate holds, an apparent 
viscosity curve measured in a parallel plate rheometer should always be smoothed out compared to one measured in 
a cone and plate rheometer. However, a comparison of measurements of Discontinuous Shear Thickening in different 
geometries showed that the appar ent viscosity curv e from a parallel plate rheometer was steeper than that measured 
with a cone and plate rheometer (|Fall et all l200st ) . This discrepancy between measuring geometries indicates that 
the constitutive relation for the stress is not a function only of the local shear rate, even for that density matched 
suspension. 

In this section we tested a constitutive relation that attributes the shape of the shear profile in the shear thick- 
ening regime for settling suspensions to the balance between viscous forces and gravitational forces for particles in 
frictional contact. The onset stress T m in occurs at the onset of particulate shear. However, the remarkable feature of 
Discontinuous Shear Thickening is the large stress jump in the shear thickening regime. We found that local consti- 
tutive relations between stress and shear rate fail to describe the stress in the shear thickening regime. Specifically, 
the collapse of the shear profile curvature values in Fig. [8] for r < T max suggests the viscous term proportional to local 
shear rate is nearly constant in the shear thickening regime, and the stress jump can not be attributed to viscous 
forces or the inhomogeneity due to gravity. Instead, the stress jump must be hidden within the additional stress term 
t c from Eqn. O whose source will be identified in Sec. [Jj 
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Figure 9: (a) A comparison of now curves measured with different boundary conditions for glass spheres 
in water at ip = 0.52 (< cj> c ). Solid circles: shear stress r for 100 /im spheres in a hxed gap measurement 
with the standard parallel plate setup. Open triangles: normal stress tn from the same measurement. The 
absolute uncertainty on the normal stress is 2 Pa, so the normal stress cannot be resolved at the low end. 
Open circles: r for 500 /im spheres with a fixed normal stress of 2040 Pa (solid line) in the modified parallel 
plate setup with a hard wall. Dashed line: slope 1 corresponding to a Newtonian scaling for reference, 
(b) Circles: same data with the shear stress r normalized by normal stress tn vs. Re — pd 2A f/r] m . Open 
triangles: constant shear rate measurements in the standard parallel plate setup in which the normal force 
was recalibrated before each measurement. Solid line: r = tjv indicating a frictional scaling. Dashed line: 
T = i)m7 corresponding to a viscous scaling. Dotted line: normal stress resolution limit. 

7 Normal forces and the boundary condition 

In this section we use measurements of shear and normal stresses under different boundary conditions to show that 
the global mechanical response can be described by a solid frictional constitutive law rather than a viscous law. 

7.1 Frictional scaling 

Here we compare steady state viscosity curves along with normal stress measurements for similar suspensions with 
different boundary conditions. The sample was 100 /im diameter glass spheres in water at a packing fraction of 
4> = 0.52 (< cj> c ). We first show results from a standard parallel plate setup (Fig. [1^,) with a diameter of 50 mm 
which results in a better normal stress resolution than smaller plates. The shear stress r and normal stress tjv are 
shown in Fig.|9K as functions of shear rate 7 for a measurement in which the gap size is fixed at 0.72 mm. The region 
with slope greater than 1 defines the shear thickening regime. We found positive normal stresses, corresponding to 
the sample pushing against the plate s, in agreement with other measurements of Discontinuous Shear Thickening 
(ILootens et all 2005; Fa ll et all l2008t ). The shear and normal stresses track each other extremely well in functional 
form and magnitude. The cutoff of tjv at the low end corresponds to the measurement dropping below the relative 
resolution of the normal stress of about 0.3 Pa. 

We next used the walled rheometer setup without a liquid air interface or room for expansion as shown in Fig. [lb. 
In this case a sample of 500 fim glass spheres in water was used; the larger particles were necessary to avoid them 
escaping through the gap between the side wall and top plate. While the values of r m i„ and T max differ with particle 
size (see Figs. f51 . otherwise the samples behave in a qualitatively similar way in the normal parallel plate setup. 
The normal force on the top plate was fixed at 1 N (2040 Pa), consequently the gap size was allowed to vary. The 
viscosity curve is shown in Fig. [9^.. In contrast to the standard parallel plate setup, the rheology is that of a yield 
stress fluid with no shear thickening regime. Such a dramatic difference in behavior with a change in boundary 
conditions would be unexpected from a local hydrodynamic constitutive relation, and implies a non-local effect. 
The common feature of both measurements is the connection of the shear stress to the normal stress. We plot the 
ratio of stresses t/tjv vs. the Reynolds number Re = pid 2A f/r] m for both measurements in Fig. |Ub. Additionally, 
we show steady state values for measurements taken at constant shear rate in the shear thickening regime in which 
the normal force was recalibrated relative to the static value before each measurement to optimize resolution of the 
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Figure 10: (a) Transient time series of shear stress r (circles) and normal stress tjv (triangles) in normal 
force control measurements for cornstarch in water at <p = 0.55 < <f> c . The sample starts at rest then the 
shear is switched on at time t — 0. Solid symbols: shear rate 7 = 26 s _1 (above j c ). Open symbols: 7 = 8 
s _1 (below 7c). (b) Same data as panel a, but extended to longer times to see the steady state behavior. 
Right axis: change in gap size Ad (small squares), (c) Effective viscosity curves obtained from transient 
measurements. Solid circles: transient shear stress averaged between 0.4 and 1.0 s after shear starts. Solid 
triangles: transient normal stress averaged over the same time. Open circles: steady state shear stress at 
the end of the time series where tn was below the resolution limit for each shear rate. Discontinuous Shear 
Thickening is suppressed when the normal stress at the boundary is removed. Solid line: stress as a function 
of shear rate obtained from a steady state viscosity curve for the same sample with fixed gap size. Dashed 
line: Newtonian scaling. 



relative normal force. The fact that these three data sets under different measurement conditions collapse onto the 
same curve suggests a global constitutive relation independent of boundary conditions. Since t/tm is near unity and 
constant over five decades of Reynolds number this suggests that the measured stresses are compressional in nature. 
A compressional scaling from either viscous, inertial, or frictional forces in dense suspensions is the result of the 
redirection of stress in different directions through the bulk of the suspension by particle interactions dNott and Bradyl . 
1994 iPrasad and Kvtomaal. 1 19951 : iBradv and Vicid. Il995l : IJaeger et all Il99d : Isierou and Bradvl. l2002l : lDeboeuf et all . 



2008T ) . This compressive scaling can also account for the constant stress term r c in the constitutive relation of Eq.(3] 



since it is independent of local shear rate and height, but instead is dependent on the normal stress. A compressional 
scaling also confirms the total normal stress is the relevant physical quantity, rather than normal stress differences. 

While viscous or inertial stresses can in some cases result in constant r/rjv, the stresses must be proportional 
to shear rate or shear rate squared, respectively, neither of which is satisfied for Discontinuous Shear Thickening 
suspensions as in Fig. [9] Rat her, a frictional exp lanation in which the forces are transmitted along chains of hard 
particles via frictional contact (I Jaeger et alj . ll99q ) is supported by the observations that the stresses have no inherent 
shear rate dependence and continue to follow this same relation when additional normal stress is applied at the 
boundary as shown in Fig. [9l the indication from the shear profile that the weight of the packing is transmitted along 
frictional contacts to build up in lower layers shown in Sec. [5J and the relative smallness of the expected viscous and 
inertial forces in the Discontinuous Shear Thickening regime, as explained in Sec. [4] 



7.2 Transient normal force control measurements 

Here we show that the coupling between shear and normal stress applies even to transient measurements as the 
normal force boundary condition cha nges, and that sh ear thickening can be eliminated if the normal force is removed 
from the boundary as suggested by iFall et al.l j2008). To emphasize the generality of these results, we show this 
result for a different suspension, cornstarch in water. Both this result and the results of the previous section were 

found for both suspensions, but we show only one set of data for brevity. 

We performed normal force controlled experiments modeled after those of lFall et al.l l|2008h . These measurements 
were done in the standard parallel plate setup in a normal force controlled mode. The normal force set point is zero 
relative to the rest state, with an initial gap of d = 1.08 mm. The gap size is free to vary during the measurements 
to adjust the normal force back to the setpoint via a feedback loop. Initially the sample of cornstarch in water at 
<j> = 0.55 < 4> c was at rest, then at time t = the shear rate was set to a constant non-zero value for the rest of 
the experiment. Examples of transient time series of the shear stress and normal stress are plotted in Fig. 110a for 
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two different shear rates. For shear rates below the onset of shear thickening -y c ~ 11 s _1 , the stress quickly came 
to near the steady state value within a fraction of a second and remained there. For shear rates above 7 C , the shear 
and normal stresses had a large peak initially, exceeding the steady state value by more than an order of magnitude. 
Even though the normal force set point was zero, the normal stress can be non-zero in the transient behavior as the 
gap adjusts via a feedback loop. Longer time series are shown in Fig. 110b along with the variation in gap size. Below 
7 C the normal force did not exceed the threshold to cause the gap to move. In contrast, above y c the gap increased 
initially due to the transient normal force. The shear stress tracked the normal stress quite well throughout the entire 
transient process, and they were similar in magnitude. The stresses each decreased as the gap increased, and the gap 
stopped increasing when the normal stress dropped below the feedback threshold of 20 Pa. Beyond this point the 
stresses and gap size remained constant, which was measured for at least 200 s in each experiment to confirm that 
the system was in a steady state. 

We summarize the normal force control experiments with effective viscosity curves in Fig. 1 10b . We show the 
transient shear and normal stresses averaged between 0.4 and 1 s after the shear was started as solid symbols. Because 
the response time of these samples to dramatic changes is typically a fraction of a second, while the normal force 
control feedback loop has a longer timescale, these transient results effectively correspond to a fixed gap boundary 
condition. They show the same qualitative shear thickening as steady state behavior for fixed gap measurements, 
indicated by the solid line. Differences between the solid circles and solid line beyond the measurement resolution 
indicate a difference between steady state and transient measurements. Stress values taken from the end of the test, 
where rjv = (±20 Pa) and the system was in a steady state, are shown as open circles in Fig. HOfc . The effective 
viscosity curve based on this data is consistent with a Newtonian scaling at shear rates below j c . Above 7 C , the shear 
stress values match up with the normal stress feedback threshold. This can be understood since a normal stress of 
that magnitude is not enough to trigger the normal force control feedback loop so the gap size remained fixed, but 
the normal stress can still couple to the shear stress. The strong shear thickening in the fixed gap and transient data 
is totally abse nt from the r j y = d ata. We note that there is no significant dependence of viscosity curves on gap size 
in this range (|Brown et alj . l2010bj ). so the difference must be due to the fact that the normal force is fixed to be zero. 
This shows that a positive normal stress of comparable magnitude is required to achieve the shear stress associated 
with shea r thick ening. In the absence of this confining stress, shear thickening cannot occur, as was suggested by 
iFall et all (120081 ). Without making any assumption about the mechanism for coupling between the normal and shear 
stresses at values below the normal stress feedback threshold, the open circles in Fig. HOfc put an upper bound on the 
viscous and other non-compressive contributions to suspension viscosity and show that they are not responsible for 
Discontinuous Shear Thickening. This directly shows the surprising result that cornstarch in water, the prototypical 
Discontinuous Shear Thickening suspension, is not actually shear thickening based on the direct constitutive relation 
between stress and shear rate when boundary conditions are held constant, and the occurrence of the phenomenon 
depends on a changing contribution of normal stresses from the boundary. 

We have noted that the normal and shear stresses track each other quite well in normal-force-controlled measure- 
ments. In fact, in all of the various types of experiments on suspensions that exhibit Discontinuous Shear Thickening 
in which normal stresses and shear stresses were compared they tended to track each other quite well. For example, 
we attempted measurements with a fixed normal stress tjv = and fixed shear stress greater than the normal stress 
feedback threshold. Since the shear stress is the dominant control parameter of the rheometer, the shear stress 
reached the set value but the measurements never reached a steady state because the normal stress c ould not drop, 



causin g the gap to increase until the top plate detached from the sample. Similarly, in experiments bv lLootens et al.l 
(|2005l ) that measured stress fluctuations in the steady state, fluctuations of the normal stress and shear stress were 
found to be strongly coupled with a proportionality close to 1. This helps explain an earlier result in which an ap- 
parent viscosity curve no longer sh owed shear thickenin g when positive fluctuations in the shear stress in the steady 
state were removed from the data (jLootens et all 120031 ) . Since the shear stress fluctuations were associated with the 
normal stress, this was in essence showing an effective viscosity cur ve with no normal str ess. Another example comes 
from our measurements of finite-size effects at very small gap sizes (|Brown et al ., 2010b), in which the normal stress 
scaled with and was close in magnitude to the shear stress as it varied with gap size. 

In this section, we showed that the shear and normal stresses are coupled with a proportionality coefficient close to 
1 in the shear thickening regime, and the relationship between stress and shear rate even changed with the boundary 
conditions to satisfy this stress coupling. These results suggest that the stresses responsible for Discontinuous Shear 
Thickening are transmitted by frictional compressive forces. 



The previous section showed that the normal stress boundary conditions controls the shear stress response. This 
implies that the normal stress must at the boundary must increase dramatically with shear under typical measurement 
conditions that produce Discontinuous Shear Thickening. The compressional scaling implies stresses are redirected 




8 Dilation 
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Figure 11: Top views of a 2.4 mm deep layer of cornstarch in water, (a): Below </> c in a shear cell at rest, 
(a) at a shear rate above j c , taken after a shear displacement of Sx = 2.5 mm relative to panel a. Dilation 
can be observed as an increase in surface roughness in the sheared region near the wall. 



through the system in all directions, so that similar confining pressures must be maintained on the system on all sides 
- even a suspension-air interface such as at the side of a parallel plate rheometer. In this section we analyze images 
of the suspension surface in contact with air to identify the boundary conditions responsible for the normal stress. 
In systems that show Discontinuous Shear Thickening, it has been long known that the r e can be a visible change in 
the surface of susp ensions at the onset of shear thickening (|Metzner and Whitlockl . 1 19581 : lO'Brien and Mackavl koOO: 
ISmith et al. . 201Ch. It w as understood early on that this was due to dilation, which occurs along with normal stresses 
( Mctz ner and Whitlockl GUI ; iRevnoldl Il885h . When wet granular packings dilate under shear, they take up more 
space than at rest, and consequently the liquid is then sucked away from the boundary into the enlarged interstices 
between grains, so by eye the surface appears to become rough as the particles poke through. 

This visible effect of dilation is shown for a suspension of cornstarch (14 (j,m) in water below <j>c in Fig. 1111 The 
suspension was in a 2.4 mm deep layer and viewed from above. One of the side walls could be displaced to shear 
the suspension. Before shear the surface of the suspension looked wet and shiny, as seen in panel a. When the 
upper right wall was sheared at a rate exceeding the onset of shear thickening the nearby suspension appeared rough, 
shown in panel b. As soon as the shear rate dropped, the surface appeared smooth and shiny again. This behavior is 
shown in supplementary video 3. We observed that the onset of visible dilatancy corresponds closely to the onset of 
shear thickening, consistent w i th previous observa tions that the shear thickening corresponds to the onset of dilation 
(|Metzner and Whitlockl. Il958l ; ISmith et alll2010l) . 

Usually we find suspensions will show Discontinuous Shear Thickening in rheological measurements if the surface 
changes from shiny to rough when sheared, indicating dilation. At low packing fractions, the surface remains smooth 
under shear because the packing fraction is too l ow for granular dilation to affect the surfac e, since volume changes 
from dilation are typically only a few percent (iRevnoldsl 1 18851 : lOnoda and Linigerl . Il990l ). Alternatively, if the 
suspension has a yield stress, the surface may be always rough and not change with shear rate, even if the packing 
still dilates with shear. Thus, the conditions where a change in the surface from dilation is observed seem to correspond 
to the conditions for suspensions to show Discontinuous Shear Thickening. 

There is a notable exception to the rule that a visible change in the surface from dilation indicates shear thickening. 
Settling p articles in a Couette ce l l wer e seen to dilate b ut did not shear thicken, and instead a yield stress was 
measured dMetzner and Whitlockl . [l95sl ; iFalT et all 12009). However, the inhomogeneity due to gravity can explain 
this behavior. The weight of the particles in a vertical column of height H pushes on the movable side wall s, resulting 
in a y ield stress on the scale of ApgH which can be on the order of kPa for typical rheometer Couette cells (|Fall et all 
120091 ) . well above the shear thickening stress regime which we observed for glass spheres from 10-100 /jm in a parallel 
plate geometry which does not measure this yield stress (see Fig. I15p . In the case of vertical walls, the side of the 
suspension is still jammed at rest which prevents shear thickening, but the top is not, so the suspension falsely appears 
unjammed when viewed from the top. A more general conclusion that applies regardless of inhomogeneities is that 
dilation is necessary for Discontinuous Shear Thickening but not sufficient because she ar thinning stresses must be 
small compared to shear thickening stresses or else the shear thickening will be hidden (|Brown et alll2010al ). 

Now that we have established the importance of the boundary conditions, we want to directly address what the 
boundary looks like on the particle scale. To this end we use a sample of opaque 150 /im diameter ZrO-2 particles 
in mineral oil at <f> = 0.54 in a standard parallel plate rheometer setup with a gap size of 890 /im. A video camera 
was focused at a point on the surface of the suspension with the line of sight tangent to the surface to best view 
any radial variations in the boundary position due to dilation. The sample is shown at rest in Fig. 112b . in which 
case it had a smooth surface. The same sample is shown in panel b at a steady shear rate of 3 s _1 , corresponding 
to T m ax at the upper bound of the shear thickening regime. The boundary appears bumpy as particles penetrate 
the liquid-air interface. The dynamic behavior is shown in supplementary video 4. By reference to the edge of the 
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Figure 12: Images a suspension of 150 /im Z1O2 particles in mineral oil in the standard parallel plate setup 
with a gap of 890 /im. The camera is focused at a point on the edge of the suspension, with the line of sight 
tangent to the surface to view radial variations in the boundary position. The rest of the image is out of 
focus because of the large amount of depth in the image, (a) The suspension at rest, (b) The suspension 
is sheared at constant shear rate of 3 Hz corresponding to T max . It can be seen that shear results in both 
radial dilation of the suspension and increased local curvature at the surface on the particle scale. Vertical 
lines: reference lines indicating the plate edge in each image. 

rheometer plate (red line), it can be seen that the sample has expanded radially relative to the rest state by about 50 
lira, corresponding to 0.3a or a volume increase of 0.8%. The penetration of the liquid-air interface by the particles 
can also be seen in the shear profile videos of Zr02 from Sec. [5] as texture differences indicating contact lines on 
particle surfaces in supplementary videos 1 and 2. 

As long as the particles are between about 1 and 100 /im, then the grains are large enough to scatter light 
diffusively and small enough that they cannot be seen individually, so the surface appears rough by eye. However, 
we note that for colloidal particles smaller than around 1 /im, the roughness of the surface becomes smaller than the 
wavelength of light and no longer scatters diffusively, so the surface of a colloid may even appear shiny if the surface 
is deformed by particles, and dilation would not be clearly visible. This is confirmed, for example, by observations 
of stable asperities in a jammed c olloid of 1.6 [im diameter part icles - indicating stresses from surface tension - in 
which the surface remained shiny (|Koos and Willenbacherl . 1201 if ) . 

In this section we suggested that the often-observed but unexplained connection between dilation and shear 
thickening is that dilation causes suspensions to interact with their boundaries and cause a change in boundary 
conditions. This could explain the observation that the boundary conditions determine the constitutive relation 
between stress and shear rate and even whether Discontinuous Shear Thickening is observed. 

9 Capillary forces 

In the previous section showed that when a dense suspension dilates under shear, the particles penetrate the liquid-air 
interface. Here we propose a model by which dilations changes the boundary condition to produce a confining stress 
from capillary forces, which can provide the normal stress required for Discontinuous Shear Thickening. Indeed, it has 
been suggest ed that capillary forces at boundaries could p lay an important role in the rheology of shear thickening 

suspensions l|Holmes et all . l200l liooi ICates et aill2005bl ). 

Ch anges in surface roughness similar to those from dilation have been observed in jammed suspensions (|Cates et all 
l2005al). and in free-surface flows, becoming more apparent at higher ap plied stresses and at higher packing fractions 
(jLoimer et all. 120021 : iTimberlake and Morrisl . [20051 : Isingh et all 12000 ). It was argued that these deformations re- 
quired normal stresses in the suspensions to balance forces from surface tension due to the curvature of the liquid-air 
interface. When dilation causes particles to penetrate the edge of the suspension to create a curved liquid-air inter- 
face, the scale of t he radius of curvatu re r of the liquid-air interface with surface tension 7 becomes comparable to 
the particle size a (|Loimer et all 12002). This produces a stress from surface tension pushing on the particles towards 
the interior of the suspension. We estimate this stress from surface tension to be on the scale of j/r ~ y/a ~ 100 
Pa for 100 /tm particles, a significant stress in the context of the rheological measurements. If the particles did not 
interact along stiff force chains, they would be pushed to the interior of the sample by this stress (assuming the liquid 
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Figure 13: The radial dilation <5 measured as a function of shear stress r for 135 fim polyethylene spheres 
in silicone oil. Solid line: predicted relationship between S and r for a model in which there is a confining 
stress from surface tension r oc j/r(6) where the local radius of curvature r(5) is calculated geometrically. A 
proportionality coefficient of 0.14 shifting the curve horizontally is used to fit the data. Dashed line: dilation 
value where the contact line is expected to reach the 2nd layer of particles resulting in a dramatic increase 
in confining stress with dilation. 



wets the particles, which is also a requirement to observe shear thickening (|Brown et al.l . lioiOal )). However, the fact 
that the particles continue to penetrate the surface in the steady state implies that, in the absence of inertial effects 
(satisfied by the low Re of the experiments), forces must be transmitted all the way through the packing to balance 
the stress from surface tension. The stress from surface tension can then be considered a confining stress which is 
transmitted and redirected through the suspension to the rheometer plates according to the compressional stress rela- 
tion seen in Figs.[9land llOI This stress would be the contribution that increases rapidly with shear rate as the system 
dilates under shear which is the characteristic feature of Discontinuous Shear Thickening. Eventually, this confining 
stress reaches the limiting scale of 7/a from surface tension, and beyond that point any additional shear stress must 
come from other sources, which are likely weak compared to the confining stress if shear thickening is observed, so 
the viscosity will likely drop off beyond the maximum confining stress. Thus the maximum confining stress should 
corres pond to T max . Abov e <f) c where the suspension is jammed, particles are seen penetrate the surface even without 
shear (|Brown et al.l . [20 111 ), so the yield stress scale Tj should also be set by the confining stress from surface tension. 
In the remainder of this section we quantitatively compare the measured shear stresses to the confining stress scale 
from surface tension 7/a to test this model. 

We now address why surface tension determines the values of r max and Tj rather than r min . We already argued 
that the confining stress from surface tension produces normal stresses between particles and the walls that results in 
a shear stress via friction. This confining stress increases precipitously in the shear thickening regime as the amount 
of dilation increases with shear rate. The maximum confining stress from surface tension should be on the order of 
7/a if 5 reaches the order of a or larger. Beyond the point where the confining stress reaches its maximum, any 
additional shear stress must come from other sources, which are likely weak compared to the confining stress if shear 
thickening is observed, so the viscosity will likely drop off beyond the maximum confining stress. Thus the maximum 
confining stress should correspo nd to Tmax- Above 4> c where the suspension is jammed, particles are seen penetrate 
the surface even without shear l|Brown et al.l . [201ll h so the yield stress scale Tj should also be set by the confining 
stress from surface tension. 

9.1 Relating dilation to confining stress 

To quantify dilation, we measured the mean radial displacement of the surface during shear seen in a tangent view as 
shown in Fig. [T2j We did this for several different steady state shear rates in a sample of opaque 135pm polyethylene 
spheres in silicone oil at tj) — 0.56. For each measurement, we started the sample at rest, then sheared at a constant 
shear rate until the stress reached a steady state for some time, then stopped the shear to observe the relaxation 
to rest. We repeated the cycle of shear followed by resting a total of 5 times. The edge of the sample was tracked 
throughout these measurements. The dilation 8 was measured as the mean radial displacement of the edge between 
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the steady state shear and rest states, averaging over the height of the sample and over at a period of at least 10 
s and a strain of at least 2 in the steady state for each experiment. The measured dilation is plotted versus the 
corresponding steady state stress values in Fig. [T3] Plotted error bars correspond to the standard deviation of 8 
measured over the 5 cycles. All points shown correspond to stresses above r m i„. At lower stresses we could not 
resolve any dilation below our resolution limit of 0.5 /im. The upper end of the shear thickening regime corresponds 
to Tmax = 2 Pa for this sample. 

Using a geometric model, we can calculate a typical radius of curvature of the liquid-air interface as it contracts 
for a given particle dilation 8, given the contact angle and conservation of liquid volume. This allows us to estimate 
a confining stress scale from surface tension r y/r(5). Details of this calculation are shown in the appendix. Briefly, 
the initial state with a relatively large radius of curvature corresponds to 5 = 0. As 8 increases, the surface becomes 
curved as particles penetrate the surface due to dilation, and the radius of curvature decreases. The corresponding 
stress scales almost linearly with 8, and the scale of the radius of curvature is set by the particle size when the dilation 
is around a particle radius. This model prediction is shown in Fig. [13] where 8 is plotted vs. the predicted stress 
scale 7/r(<5). A free parameter for the scale factor of 0.14 on the stress scale is used to fit the data. The qualitative 
agreement in the model slope with the data in Fig. 1131 confirms that the confining stress scaling as ■y/r is a good 
estimate for the measured shear stress. The fit coefficient within an order of magnitude of 1 confirms that the amount 
of dilation is on the right scale to provide the measured stress. 

The dotted line in Fig. [13] corresponds to the dilation value where the contact line is expected to reach the 2nd 
layer of particles from the surface (see appendix for calculation). At this point the confining stress should increase 
rapidly as more contacts are made with small curvature. Because the calculation of confining stress from dilation is 
not monotonic around this region, the dilation is not single-valued function of confining stress. The lower portion 
of the curve is expected to be unstable since more dilation would provide less of the stress required to confine the 
suspension to a smaller volume. The agreement of the dilation measurements with the dotted line beyond the point 
where the calculation becomes multi- valued supports this interpretation. 

As the dilation increases and the contact line recedes further into the interior, a lower limit for the value for the 
curvature must be reached as it is limited by the interstitial gap size. The corresponding limiting confining pressure 
has been measured in an analogous system in which a fluid interface was driven through a porous m edium, in which 
case the required driving pressure went to 0.7j/a in the limit of zero flow rate (|Weitz et all Il987h . This confining 
stress sets the scale for the upper bound on the component of the shear stress due to capillary forces on the order of 
y/a in the limit of large 8, although the exact value of the coefficient should depend on material specific properties 
such as contact angle and particle roughness. The value of t mi = 2 Pa corresponding to the data in Fig. [13] is 
significantly below the limiting confining stress regime, suggesting that dilation by a fraction of a particle width was 
enough to obtain a fully developed shear flow and the limiting confining stress is not necessarily reached in the shear 
thickening regime. 

9.2 Surface tension scaling 

To confirm the role of capillary forces, we performed a set of rheological measurements in which we varied the surface 
tension of the liquid-air interface. To vary this surface tension, we added surfactant to 100 fim glass spheres in 
water. The surfactant used was Palmolive dish detergent, which was first mixed in water above the critical micelle 
concentration which reduces the surface tension with air by about a factor of 3 compared to pure water and air. 
Viscosity curves are shown with and without surfactant and at different packing fractions in Fig. 1141 We first 
compare the viscosity curves for jammed suspensions at <j> — 0.58 > 4> c . These viscosity curves correspond to yield 
stress fluids. The value of the yield stress is reduced by a factor of 2.4 with the addition of the surfactant, about the 
same as the surface tension was reduced. 

We next compare the viscosity curves at <j> = 0.56 < <f> c in Fig. 1141 The increase in the viscosity at low shear rates 
can be attributed to the increase in the zero shear viscosity with the addition of the surfactant. In terms of stress 
scales, there is a decrease in T max by a factor of 2.4 when the surfactant is added, and no resolvable change in r m ; n . 
The reduction in both Tmax and the yield stress r y above <j> c is comparable to the reduction in surface tension with 
the addition of surfactant, again consistent with a model in which these stresses scale with surface tension. 

We note that in principle the addition of surfactant can change other relevant parameters. The stress from surface 
tension on a boundary typically scales as (7/r) cos# where r is the radius of curvature and 8 is the contact angle, 
where both 7 and 9 can vary with the addition of surfactant. The addition of a surfactant can reduce 8, increasing 
the stress from surface te nsion. However, we st art with a liquid that wets glass pretty well, as this is a requirement to 
observe shear thickening l|Brown et al. , 2010a), so cos 9 ~ 1 even before the addition of surfactant. The addition of 
surfac tant can also affect the value of r m i n in cases where the particle- liquid surface tension is dominant (|Brown et all 
l2010al 1. but for these 100 /im glass spheres in a wetting liquid the dominant force affecting the onset stress is gravity 
(Fig.©. 

The results in Fig. [14] suggest that the upper stress scales T max and Tj (the scale of the yield stress for <j> > <j) c ) 
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x(Pa) 

Figure 14: Viscosity curves for 100 /Ltm glass spheres in liquids with different values of surface tension. 
Solid symbols: particles were suspended in water. Open symbols: particles were suspended in water with 
surfactant (above the critical micelle concentration). Triangles: (f> — 0.58 > <p c . Circles: 4> = 0.56 < </> c . 
Solid line: r max without surfactant. Dotted line: T max with surfactant. Both r max and the yield stress above 
cj) c decreased when the surface tension was reduced. 
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Figure 15: The stress at the upper bound of the shear thickening regime T max for a variety of suspen- 
sions plotted against the confining stress scale from surface tension 7/a. Particle mate rials are l i sted in 
the key. Solid symbols: measured by us. Open symbols : polyvinyl chlo ride [PVC, circles ([Hoffman . 1972 )]. 
polystyrene-acrylonite [PSAN, down-p ointing triangles (iHoffmanl . I1972I) ]. polys tyrene [up-pointing triangl es 
( Boersma et al. . 199 lh ]. g lass [square dBoersma et aL , 1990 )]. silica [diamond dBender and Wagner . 1996 )]. 
CaCO s [diagonal crosses (lEgres and WagnerL 120051 )]. PMMA [cr ossed square jKalman et all 2009)]. BiOCI 
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scale with the surface tension at the liquid-air interface. To more generally test the predicted stress scale T max ~ 7/0 
including the particle size scaling, we plot measured values of T max vs. 7/a for the wide range of Discontinuous Shear 
Thickening suspensions we have studied in Fig. 1151 Each point corresponds to a different suspension, with a wide 
range of different particle materials, shapes and sizes, and different liquids. The plotted value of T max is an average 
over viscosity curves at several packing fractions. We also included data from other papers in cases where T ma x was 
measured. It is seen that for this wide variety of suspensions, covering four orders of magnitude, r mal falls in a band 
that scales as the prediction 7/a indicated by the solid line. We note that for each Discontinuous Shear Thickening 
suspension we studied, the two stress scales Tmax and t 7 - are alw ays within an order-of-magnitude of each other, as 
was seen, for example, in Fig. 1141 and Bro wn and Jaeger] (|2009l ). suggesting that Tj also scales with 7/a. In many 
measurements of colloids, the upper end of the shear thickening regime was not reached. If there is an upper bound, 
it would have to be above the range measured. This is especially a problem with colloids because the expected scale 
of Tmax for small particles exceeds the measuring range of many rheometers. For example, our Anton Paar MCR 301 
rheometer has an upper limit of 3800 Pa for the Couette cell or 65,000 Pa for the 25 mm diameter parallel plate. 
This lowe r bound on T max based on the l imited measuring range is illustrated as the dotted line in Fig. 1141 using 
data from M aranzano and Wagner! (|2001aT ) as an example. 

There is variation in the value of T ma x in the band shown in Fig. [15] by about an order of magnitude. There are 
many factors that could contribute to the precise value of the confining stress and the resulting shear stress. For 
example, the normal stresses do not have to be exactly the same on each surf ace as would be the case fo r a pressure 
acting on a fluid. Instead the stresses are related by a coefficient of order 1 IjJansse 3, 1 18951 ; ISperll. |2006ft . Since the 
confining stress can put a normal stress on the top plate via chains of particle contacts, then a component of the shear 
stress comes from friction, related to the normal stress by an effective coefficient of friction as seen in Fig. [9]which can 
depend on many factors including particle shape, roughness, and particle interactions. The contact angle 9 has been 
left out of the force equation since it is not known in many cases. The dependence of dilation on the shear rate must 
also play a significant role, as seen in Fig. 1131 Geometric factors including particle shape and roughness also should 
play a role that has not yet been studied. Considering all of these dimensionless factors of order 1 that can affect the 
shear stress which are not all known or easily measured, we will not go beyond using the order-of-magnitude stress 
scale of 7/a as an estimate for T max . 

In this section we proposed that Discontinuous Shear Thickening is due to a confining stress from surface tension 
in response to deformation of the liquid-air interface at the boundary from dilation. We confirmed the stress increases 
with dilation, the effect of surface tension, and showed that r mal scales with the confining stress scale 7/a from surface 
tension with more than 30 suspensions that cover four decades of particle size. 

10 Solid boundaries 

In the previous section we showed that under boundary conditions such that particles penetrate the liquid-air interface, 
surface tension provides the confining stress that is responsible for Discontinuous Shear Thickening. While a liquid-air 
interface at the boundary is typical for rheometer measurements, closed systems with solid walls are also of interest. 
In this section we will use a solid-walled rheometer setup to investigate the role of the confining stresses in closed 
systems. 

For measurements in a closed system, we used the parallel plate setup with solid walls shown in Fig. [TJs. The hard 
walls confined large grains within the container volume without the need for the surface tension of the liquid. Thus, 
we can also determine the role of the liquid by comparing measurements with and without liquid. We first show 
stress vs. shear-rate curves for shear rate controlled measurements of dry 500 /im glass spheres in Fig. 1161 Without 
liquid, the packing fraction is determined by the container volume which can be varied with the gap size. Thus, for a 
series of measurements with a fixed volume of particles, the gap height determines the packing fraction, with smaller 
gaps corresponding to higher packing fractions. We give packing fraction values accurate to 3 decimal places relative 
to each other to compare curves in Fig. 1161 but absolute uncertainties on packing fractions are still around 0.01. 

An important result from Fig. [16] is that the curves show Discontinuous Shear Thickening that is qualitatively 
similar to measurements of suspensions in standard rheology setups, despite the fact that there is no liquid. Thus, the 
interstitial liquid or viscous stresses are not a necessary component for shear thickening when the grains are confined 
by other means. 

A large hysteresis loop can be seen for <f> = 0.625 in Fig. 1161 This is the threshold beyond which - at larger gap 
sizes, corresponding to lower packing fractions - not only the yield stress but also the measured shear stress was 
below the resolution limit suggesting contact between the plate and grains was lost. This emphasizes that a key role 
of the liquid is simply to keep contact with the plates and transmit stress between the particles and the plate. 

We repeated these measurements with water as a solvent filling the measurement volume and the surrounding 
volume so there was no liquid-air interface near any particles. With water, contact between the suspension and 
plates could be maintained at larger gaps (lower packing fractions). However, no significant difference was seen in 
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Figure 16: Stress vs. shear rate curves for 500 /im diameter glass spheres in a solid-walled rhoemeter with 
no liquid. Packing fractions <j> shown in the key; higher curves correspond to larger <p. Discontinuous Shear 
Thickening is still seen, confirming that viscous interactions are not necessary. 

the qualitative aspects of Discontinuous Shear Thickening or in the scale of T max with or without water. Notably, 
the scale of r max with the solid wall is almost two orders of magnitude higher than for the same suspension in the 
parallel plate setup with a liquid-air interface (Fig. I15|l . 

An upper bound on the jamming transition can be identified by the point where the shear stress drops below the 
measurement resolution at <f> = 0.62. The jamming transition can be at a significantly higher pack ing fraction dry than 
with liquid because of the larger density difference between the particles and surrounding fluid (|Onoda and Linigerl . 
Il990h . With a solid wall, the yield stress did not plateau at high packing fractions like in the case of a liquid-air 
interface, but rather increased dramatically as the packing fraction was increased as seen in Fig. 1161 This continued up 
to the maximum stress the rheometer can apply. This can be expected if the confining stress comes from the stiffness 
of either the wall or the particles, in which case the confining stress increases as the solids are further compressed 
(|Q'Hern et aUl2003h . This is in contrast to the confining stress from surface tension w hich is limited by the cap illary 
stress which reaches a plateau value as packing fraction is increased beyond jamming l|Brown and Jaeger] , 120091 ) due 
to the fact that the minimum radius of curvature is set by the particle size. 

To connect Tmax to a confining stress for closed systems, we made measurements of the compressional stiffness of 
the tool and sample in series. We observed that the rearrangement of particles under shear makes the suspensions 
much more compliant than under compression alone. Thus we sheared the samples while measuring the compressional 
stiffness to better match the usual experimental conditions. The sample was slowly compressing at a fixed rate of 
0.25 /im/s while also shearing at a fixed rate of 1 Hz (3 mm/s). The shear rate was much faster than the compression 
rate so that the packing has time to rearrange as it is being compressed, but slow enough that the shear stress is 
still near the zero shear rate limit as seen in Fig. 1161 The measured shear and normal stresses are shown in Fig. [17] 
We note that the stiffness under shear is much less than the value obtained by compressing the sandpaper by itself, 
which is the weakest component of the wall. Thus the presence of the grains has a significant effect on the effective 
stiffness, despite the fact that the material stiffness is much higher than that of the sandpaper. 

In the case with a liquid-air interface, we found the confining stress to scale roughly as —S'y/a 2 due to linear 
compression of the boundary from dilation (Fig. [15) . If for the solid walls, the confining stress also comes from 
compression of the boundary, it should provide a restoring stress in response to dilation of Sdr/dd. To obtain an 
analog for the surface tension so that the confining stress scale can be written as 5k/ a 2 , we define a stiffness per 
particle as k — —a 2 dr/dd. This differs from the usual definition of stiffness for an elastic material; rather than being 
proportional to wall surface area, it is normalized for a wall whose cross-sectional area is a 2 , near that of a particle. 
Since this normalization for the per-particle-stiffness makes it independent of the system size, it is more physically 
relevant in discussions of stress scales. We can obtain k from the slope in Fig. 1171 

We plot values of T max vs. a confining stress scale equal to the boundary stiffness per particle over particle size 
k/a in Fig. 1181 each for the wet and dry 500 (im glass spheres. We also measured a set of stress/shear-rate curves 
and compression curves for the wet glass spheres with a layer of soft foam rubber inserted between the top plate 
and sample as shown in Fig. []J>, and the corresponding values of T ma x and k/a are also plotted in Fig. [18] We also 
measured a set of data without sandpaper on the plate surfaces, which resulted in a much lower stiffness than with the 
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d (mm) 

Figure 17: Shear stress t (open triangles) and normal stress tn (solid circles)vs. gap size d for a sample of 
500 ^m glass spheres with no liquid under slow compression with solid walls. The sample is compressed at 
a rate of 0.25/im/s and sheared at a rate of 1 Hz (3mm/s). The shear stress is close to the zero shear rate 
limit, so the measured r is a good proxy for the yield stress. The solid line is a linear fit used to obtain the 
per particle stiffness k of the system of sheared grains and solid wall in series. 
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Figure 18: Maximum stress of the shear thickening regime T max vs. the confining stress scale k/a due to the 
restoring force of the boundary with per-particle stiffness k. Data are for 500 glass spheres under several 
different boundary conditions. Solid triangle: hard wall rheometer setup, with particles suspended in water. 
Open triangle: hard wall, no liquid. Solid square: hard wall with a soft foam rubber insert, with particles 
suspended in water. Open circle: in a standard parallel plate setup with a liquid-air interface, where we use 
the surface tension to represent the per particle stiffness (i.e. k = 7). Solid circle: polyethylene in mineral 
oil from Fig. [T3l where we calculate k — a 2 (85 /dr)^ 1 . The solid line has a slope of 1, corresponding to a 
stress response proportional to the restoring force of the boundary against a typical dilation of the sample 
by 6 w 0.05a. 
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sandpaper, despite having a harder surface, suggesting significant slip. For comparison to the standard parallel plate 
measurements with a liquid-air interface, we plot T max vs. y/a for one such experiment with the same particles, using 
the surface tension 7 as a proxy for the stiffness per particle, which was also plotted in Fig.[T3] As a test comparing to 
directly observed dilation, we plot T max vs. an effective stiffness per particle taken from Fig. I13l as k — a 2 {d8 / Ot)^ 1 
in the nearly linear regime. All of these experiments are consistent with the relationship Tmax = 0.05fc/a. This scaling 
confirms that, for a wide range of boundary conditions including both liquid and solid boundaries, and even when 
there is a large amount of slip, the scale of Tmax is set by the confining stress which is proportional to the per particle 
stiffness of the boundary. The similar values obtained when using k/a or 7/a as the effective stiffness confirms that 
the per-particle stiffness k generalizes the role of the surface tension with a liquid-air interface to the case with a solid 
wall. This means that the slopes of Figs. [15] and [18] are related, although with the value of the coefficient relating r ma x 
and k/a being suspension dependent. The comparison with direct measurements of dilation from Fig. 1131 confirms 
that the confining stress could be written directly as r con f = 8k j a 2 for a linear elastic boundary. 

Since k/a represents the restoring stress from the boundary for <5 = a, and the coupling coefficient between shear 
and normal stresses tends to be close to 1 (Fig. 0, the slope of 0.05 between T max and k/a suggests the restoring 
stress of the boundary is against a dilation of approximately 8 ~ 0.05a at T max . This is the order of the expected 
dilation required to mobilize a single layer of shear, for example to allow spherical particles to escape out of the wells 
they can sit it at the interstices between three neighboring particles. Similarly, dilation on the order of a percent of 
the s ample thickness is typical of measurements of sheared granular packings l|Revnoldsl . Il885l ; lOnoda and Linigerl . 
Il99d ). 

In this section, we showed that for closed systems, the maximum stress in the shear thickening regime T max is 
determined by the wall stiffness. In the experiments with a liquid-air interface, the stiffness of the boundary at the 
side determined the confining stress, while in the solid wall experiments, the stiffness of the top wall with the soft 
layer determined the confining stress. In each case, the stiffness of the most compliant boundary determined the 
confining stress. This is due to the fact that the stiffness of a collection of materials in series is determined by the 
most compliant material. Whether the most compliant boundary is on the side or the top does not matter, since 
the compre ssional relation between shear and normal stresses suggests that the stresses are similar on all of the 
boundaries l|janssenl . [l895l ). 

11 Discussion 

11.1 Non-local constitutive relation 

In this subsection, we combine the results of the previous sections to write the full constitutive relation for the 
mechanism we have found for Discontinuous Shear Thickening suspensions. In Sec. [6] we showed a constitutive 
relation for the shear stress accounting for viscous and gravitational forces. The shear profile measurements implied 
that the dramatic increase in stress associated with Discontinuous Shear Thickening could not be attributed to a 
local shear rate dependence from viscous stresses, whose contribution appeared to be relatively constant in the shear 
thickening regime. Instead, we showed in Sec. [7] that the shear stress was coupled to the normal stress, with a 
proportionality given by an effective friction coefficient /x. In turn, the normal stress was shown to come from the 
restoring force of the boundary in response to dilation (Sees. [9] I10| l. We can now fully express the constitutive 
relation for local shear stress tj based on Eq.[3]with the uniform term r c explicitly separated into the confining stress 
contribution as well as a constant contribution from other interparticle interactions Tint'. 

n = Tj„((j>)jif(Ite) + Tgh/d + [iT conf (8) + r int ■ (8) 

The stress terms represent, from left to right,: viscous and inertial hydrodynamics, gravity, confinement, and particle 
attractions. The confinement term simplifies to T con f(8) = 8k/ a 2 for a linear elastic boundary with per-particle 
stiffness k in response to dilation 8 (Fig. I18[) . The hydrodynamic term varies with Reynolds number such that 
/(i?e) = 1 in the viscous regime where Re < Re c ~ 100 and f(Re) oc Re/Re c in the fully inertial regime where 
i?e> Re c ~ 10 2 (Sec. U. 

While this constitutive relation is approximate, its main use is that it allows for simple estimates of the stress 
scales T m i n , Tmax , and Tj in different regimes, and their comparison can easily determine whether a system will shear 
thin or thicken. This includes, for example, the requirement that the increase in stress in the confine ment term in 
respon se to dilation must be large compared to any yield stress from the attraction or gravity terms (IBrown et all 
l2010al ). It also delineates the packing fraction regimes that separate inertial and Discontinuous Shear Thickening as 
seen in Fig. [3] The difference between the boundary conditions of fixed gap and fixed normal force as in Fig. [9] can 
be characterized by whether the confinement term is fixed by the boundary or varies in response to dilation. 

The local constitutive relation can also simply explain observations under different measuring geometries, for 
example in suspensions where gravity is dominant. In a parallel plate geometry, we found settling suspensions to 
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exhibit an apparently Newtonian scaling behavior below r m i„ because the plate measures the shear stress at a depth 
h = 0. In a Couette geometry wh ere the moving w all is on the side, any density mismatch was found to result in a 
yield stress scaling as r g ~ ApgH (|Fall et all |2009| ). since the average depth at the moving wall is (h) = H/2. 

It is interesting to note that the direct shear-rate dependence in Eq.[8]is inherently shear thinning in the low- .Re 
regime where Discontinuous Shear Thickening is observed. Similarly, we found that the local viscosity in the bulk 
obtained from shear profile measurements corresponded to shear thinning (Figs. I6l7)) . even for a global mechanical 
response that corresponded to shear thickening. The explanation for this apparent contradiction is that most of the 
shear stress is due to frictional contacts and normal stress in response to frustrated dilation, which only indirectly 
depends on shear rate. This response is non-local in the sense that the response depends on the global dilation and 
boundary conditions. Dilation can be treated as a global rather than local parameter as long as the boundary is 
linearly elastic, which is a good approximation for most solid boundaries and for surface tension where k — 7 (Fig. 1131) . 
so the average normal stress from the boundary is proportional to the average strain from dilation, regardless of 
inhomogeneities. Local variations only need to be accounted for if the boundary stiffness is non-linear. 

It is notable that the boundary conditions can have an uncommonly large effect on the normal stress term 
because it is due to a confining stress. From a hydrodynamic point of view, the large significance of the boundary 
conditions and difference between local and global results is unusual, as it is more typical for the boundary conditions 
to play a smaller role, requiring only perturbative corrections to translate between the local and global rheology. 
One of the surprising consequences of this is that characterizing rheology based solely on local, shear-rate dependent 
constitutive laws or local viscosities in the bulk would miss the dramatic phenomenon associated with Discontinuous 
Shear Thickening. 

We cannot yet fully solve the constitutive relation in Eq. [8] because dilation in suspensions is not yet well- 
characterized in the Discontinuous Shear Thickening regime of concentrated suspensions. Specifically, the slope of 
viscosity curves in the shear thickening regime should depend on how dilation couples with shear rate, stresses, and 
packing fraction. Dilation has been char acterized in suspensions tha t are dominated by viscous interactions and ar e 
not so confined as to frustrate dilation l|Prasad and Kvtomaal . fl995l ; ISierou and Bradvl . |20o3 : iDeboeuf et at [2008). 
On the other hand, d ilation in dry grains does not occur until significantly higher packing fractions, at least in 
the quasi-static limit |Onoda and Linigerl . Il990l ; I Jerkins et all 120081 ; iKabla and Sender! I2009T ) . We speculate that 
the hydrocluster models may be able to bridge this gap by describing the transition from viscous flow to frustrated 
dilation with fricti o nal c ontacts l|Q'Brien and Mackavl . l200Ch . Perhaps if hydrodynamically-induced particle clusters 
(|Bradv and Bossisl . Il985h were to become large enough to span the system become jammed against the boundaries, 
they could lead to the frictional contacts and dilation. 

The shear profile measurements of Sec. [5] suggests a possible connection between shear banding and dilation. 
Fig. 1 131 shows a single scaling for the dilation with shear stress through the transition at r max . Since the stress/shear- 
rate curve corresponds to Discontinuous Shear thickening, this implies the dilation 8 increases rapidly with shear rate 
in the shear thickening regime, then increases less rapidly above T ma x- Suggestively, this rapid increase of dilation 
with shear rate corresponds to the regime where the shear band is widening (Fig. [6]), and the dilation increases less 
rapidly with shear rate once the shear band stops widening. We cannot say for sure that these observations are 
connected, but it may be that the rapid increase in dilation is required to involve more layers of particles in the shear 
flow. 

11.2 Dominant stress scales for the onset 

In this subsection we discuss the different scaling laws for the onset stress T m in found in different parameter regimes, 
including suspensions dominated by gravity and colloids dominated by Brownian motion and electrostatic forces. We 
suggest that different scalings for r m i„ can be understood as the dominant stress scales of particle interactions in 
each system, and these observations suggest a more general requirement for the onset of shear thickening. 

We first compare the different scaling laws for the onset stress r m i„ to identify common features of the different 
regimes. In Sec.0 we found that for suspensions of particles large enough to settle the scale of T m in is set by gravity 
such that the onset of shear thickening required enough shear stress to lift the weight of the top layer of particles 
to initiate shear and dilation (Fig. [5]). In a prev ious work we conside red the effect of induced particle attractions in 
response to external electric and magnetic fields (|Brown et all [2010al V In each case the attractions resulted in a yield 
stress. The scale of r m i„ was set by the shear stress required to overcome roughly the two-particle attractive force 
(per cross-sectional area of a particle) to shear them apart. In both the gravity and attraction regimes the shear stress 
must simply exceed all local stress barriers from any source that are responsible for preventing relative shear between 
particles. While this generalization is sufficient for frictional systems, hydrodynamic systems can exhibit shear at 
stresses significa ntly below the onset of shear thi ckening. One such regime is colloids stabilized by an electrostatic 
zeta potential f. iMaranzano and Wagner! |2001al ) measured T m in for such particles from 80 to 700 nm in diameter. 
A power law fit to their measurements for the onset of shear thickening gave T m i„ oc a _211±0 ' 16 . An electrostatic 
calculation of the two-particle repulsive force over the cross-sectional area of spherical particles gives a stress scale of 
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16e£ 2 /a 2 for a liquid permittivity e. This has both the same scaling and magnitude within about a factor of 2 of their 
data, consistent with the idea that the scal ing for th e onset of shear thickening is determined by the interaction stress 
scale. The same scaling was argued for by Hoffman (1998), although with a somewhat different mechanism in mind. 
In another regime dominated by B rownian motion, the onset of shear thickening has been found to c orrespond to 
an onset stress r m i n — 50kT/3ira 3 l|Maranzano and Wagnerl . l2001bl ; iGopalikrishnan and Zukoskil [2004h . This scale 
is the osmotic pressure, which is effectively the pressure which neighboring particles interact with. The common 
feature In each of these scaling regimes is that the onset of shear thickening corresponds to the dominant stress scale 
at which neighboring particles interact, whether they interact by induced dipoles, gravity, zeta potential, or osmotic 
pressure, in each case. With particular relevance to mechanisms for shear thickening, these stress barriers prevent 
compressive shear of particles into each other, and once these stress barriers are overcome, compressive interactions 
are possible which could lead to dilation against a confining boundary to produce Discontinuous Shear Thickening. 

It is notable that the scaling laws for the onset stress r m j n have been found to work both for Discontinuous 
and Continuous Shear Thickening (Ma ranzano and Wagnerl. l2001bl ; IGopalikrishnan and Zukoskil . 120041 ). even though 
in some cases the mechanisms may be different. We suggest the generality of the onset law arises because the 
scalings for T m in come from interaction scales that are not directly responsible for shear thickening, rather they are 
preventing it. This is why different proposed models h ave predicted the same scaling la w for r m i n in the zeta-potential 
domina ted regime, incl uding the hydrocluster model (|Maranzano and Wagnerl . l2001aT ). the order-disorder transition 
model ([Hoffman! . 1 19981 ). and the dilational model, even though the proposed mechanisms for shear thickening are 
very different. The hydrocluster model which describes Continuous Shear Thickening also requires compressional flow 
between particles to occur, and thus the same stress scales would prevent both mechanisms. Thus, the requirement for 
shear thickening that the shear stress exceed the various stress scales of particle interactions that prevent compressive 
shear seems to be valid for both Continuous and Discontinuous Shear Thickening, in both suspensions and colloids, 
and whether the particle interactions are frictional or hydrodynamic in nature, as long as the mechanism for shear 
thickening requires compressional interactions. 

Many of these sa me stress scales that determine the onset of shear thickening can also determine the yield stress 
(|Brown et alj . l2010al ). with the exception of osmotic pressure. Just above the yield stress, the confining stress can start 
to grow as shear causes dilation, but there must b e at least a small shear thinn i ng regime before the confin ing stress 
becomes dominant and shear thickening is seen (|Gopalikrishnan and Zukoskil |2004| ; iBrown et al. I. l2010al ). Indeed, 
iMetzner and Whitlockl (| 19581 ) found the onset of dilatancy at slightly lower shear rates than the onset of shear 
thickening. In other words, the onset of shear thickening corresponds to a transition in the dominance of different 
stress and not necessarily where the mechanism for shear thickening first appears. 

While we have described a mechanism for Discontinuous Shear Thickening that is based on generic phenomena 
such as dilation, not all suspensions and colloids exhibit this behavior at high packing fractions. This can be explained 
by the relative importance of different stress scales. If any other particle interaction scales exceed the c onfining stress 
from s urface tension, we would expect shear thinning mechanisms to be dominant over shear thickeni ng (IBrown et al. . 
2010a). Some of the interactions that we did not exp licitly mention above include hydroge n bondi ng (Ragh avan et al. . 
2000), depletion (IGopalikrishnan and Zukoskil . 120041 ) . or a particle- liquid surface tension (|Barnesl . ll989l ; lBrown et aL . 
2010al ). One reason that many dense suspensions and colloids do not exhibit shear thickening is that many of these 
systems fall into regimes where other stress scales exceed the confining stress scale so they do not have any observable 
shear thickening regime. 

Another example of shear thinning stresses hiding shear thickening can be seen in Fig. [S] for settling suspensions. 
The maximum particle size at which shear thickening was found was about 1000 /jm. At such a large scale, the 
weight of the particles contributes to a large frictional stress r m i„ ~ Apga, which increases with particle size. This 
can overwhelm the scale of the confining stress from surface tension T max ~ j /a, which decreases with particle 
size. The balance between these two stress scales occurs at a particle size a ~ -y/7 /(Apg) ~ 1000 jj,m, above which 
the shear thinning stresses dominate, in agreement with the maximum size particle found to shear thicken. This 
transition scale could be interpreted as a particle capillary length scale which differs from the usual capillary length 
in two ways. First, this particle capillary length depends on the density difference rather than just a liquid density, 
so density matching could allow for shear thickening of larger particles. Second, this particle capillary length depends 
on particle size rather than system size. This means surface tension effects can be seen in suspensions on much 
la rger scales than wou ld be expected based on the usual capillary length. A similar result was found by the work 
of|L oimer et all l|2002h on free-surface flows of dense suspensions, in which they found effective stresses from surface 
tension scaling as "//a. 

11.3 Connection to confining stresses in other systems 

In this subsection, we discuss similarities of the role of the confining stress to other systems, including colloids, soil 
mechanics, and jammed systems. Similar to the case for T m in, we also suggest that different scalings for T max can be 
understood in a more general framework. 
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We found the upper end of the shear thickening regime r max to be set by the maximum confining stress which 
comes from the restoring force when grains dilate against a boundary, either from surface tension when there is 
a liquid-air interface (Figs. 1141 I15|) or by the stiffness of the wall when all boundaries are solid (Fig. I18p . More 
generally, this role of the confining stress requires only compressional interactions between particles that lead to 
dilation, whether the forces are transmitted via frictional interactions as we showed in Sec. [JJ or through lubrication 
interactions as suggested by hydrodynamic models. The importance of the confining stress can also be independent 
of what mechanism initiates shear thickening as long as the particle interactions are compressive. 

In the discussions of the confining stress so far, the most compliant boundary set the response. This is because 
the stiffness of a system of several elements with very different stiffnesses in series will generally be determined by 
the most compliant element in the series. In some systems, the particles could be the most compliant element. 
This regime would be relevant when all of the boundaries are hard in comparison to the particles, and it has been 
proposed such a regime may be rea ched for small colloidal particles where the confining stress from surface tension 
is larger (Wa gner and Bradvl . 2009). If there is a lubrication layer of liquid between gaps, the maximum confining 
stress would be coupled to the viscosity because the particle compression depends on the stress in the lubrication 
layer. For solid contacts between elastic spheres, the confining stress would be limited by a scale of (5/ R) 3 ^ 2 E p where 
E p is the compressional modulus of the particles and S/R corresponds to the compressional strain on the sample. 
The 3/2 power comes from the contact between two spherical surfaces as opposed to the power of 1 for flat surfaces. 
For the hard particles we used with E p ~ 10 10 and 8 ~ 10 -2 , this scale is of order 10 7 Pa, which is much stiffer 
than the liquid-air interface. Stresses up to about 1 7 Pa have been observed in the shear thickening regime for 
silica particles in compressional flows with solid walls dLim et all |2010| ). which is consistent with the idea that much 
larger confining stresses can be reached for hard walls and particle stiffness becomes the limiting factor. On the other 
hand, if the particles are extremely soft, the confining stress could be below r m i„. This seems to be the case in a 
suspension of soft gel particles with modulus on the order 10 4 Pa in which only shear thi nning was observed instead 
of shear thickening as the jamming transition was approached l|Nordstrom et all . |2010| ). So far, it is not certain 
whether or not Discontinuous Shear Thickening is controlled by a confining stress in the colloid regime. Despite the 
comparison between measurements of T max and the confining stress from surface tension in Fig. [TS] and the arguments 
for a particle-stiffness limited confining stress, there is still no direct experimental evidence with a systematic control 
of boundary or particle stiffnesses to conclusively determine whether Discontinuous Shear Thickening in colloids is 
limited by a confining stress. 

While simulations of suspension have successfully modeled Continuous Shear Thickening effects, so far most 
have failed to produce t he large, steady-sta t e stress in c reases associated with Discon t inuous Shear Thickening 
(iBradv and Bossisj . Il985l ; IMelrose et ail . Il996l ; iFarr et al l . 1 19971 ; iBergenholtz et al l . |2002| ; iMelrose and Ball l2001al . 
l2004bl ; Greben kov. 2008). These simulations have included viscous interactions as well as various interparticle in- 
teractions. Most have focused on bulk behavior, usually using periodic boundary conditions such as Lees-Edwards 
to avoid dealing with boundary effects. Now that we have recognized that the boundary conditions and especially 
the confining stress are important, it seems likely that many of these simulations did not find Discontinuous Shear 
Thickening because of their treatment of the boundary conditions. The one simulation we are aware of that has 
produced Discontinuous Shear Thickening was a molecular dynamics simulation of t wo-dimensional granular shea r 
flow with frictional contacts between particles but no liquid or viscous interactions (lOtsuki and Havakawal . [2010). 
Besides a steep T(j) at packing fractions just below the jamming transition, the scale of the normal and shear stresses 
was found to be set by the particle modulus, which in that simulation was the only scale that could limit a confin- 
ing stress. This may be a minimal model for Discontinuous Shear Thickening in two dimensions since it includes 
particle-particle contacts with a restoring force, but leaves out the liquid. 

While we have described a mechanism for shear thickening due to confining stresses in shear flows, the same 
principle could apply to extensional and compressional flows because stresses tend to be easily redistributed in 
different directions when the shear and normal stres ses are proportion al. Visible dilation at the surface was seen to 
correspond to shear thickening in an extensional flowdSmith et alll201Ch. The maximum stress in the shear thickening 
regime T max in extensional flows (jChellamuthu et aU 20091 ; Bischoff White et all l2010h has been found to be about 
an order of magnitude higher than T max based on shear measurements, but this could still be consistent with a scale 
of 7/a. While these isolated result are promising, our model has not yet been extensively tested in extensional flows. 

Systems where the mechanics a re determined by a confining stress are already well-known, for example, in granular 
systems, especially soil mechanics l|Lambe and Whitman! , [l 969), One lesson to take away from soil mechanics is that 
even though the global response of the system is set by the boundary conditions, the scale of the stress response is not 
dependent on sample size or shape. This is because forces will transmit throughout the bulk across particle contacts, 
and forces must balance across the system, regardless of how far across the bulk is. This makes stress the appropriate 
size-independent force scale as in other continuum systems. In Discont inuous Shear Thicke ning suspensions, the 
same qua litative behavior has been seen from as few as 2 particle layers l|Brown et all l2010bl ) to tens of thousands 
of layers l|Maranzano and WagnerllioOlal ). Quantitatively, the significance of the surface area to volume ratio can be 
checked by varying the gap size in a parallel plate geometry for a fixed volume of sample. For such measurements at 
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constant shear rate, w e found that the perc entage change in stress was (0.20 ±0.22)% over a range where the surface 
area changed by 17% (|Brown et all l2010bl ). This is consistent with a shear stress independent of surface area and 
inconsistent with a stress proportional to surface area. Furthe rmore, similar re sults have been found when comparing 
Couette cell measurements with parallel plate measurements (|Fall et aU . l200a ). Together, these observations support 

the argument that the relevant stress scale for shear thickening is not dependent on t he system s i ze or shape. 

It has been suggested that Discontinuous Shear Thickening is a form of jammin g l|Farr et alj . 1 19971 : ICates et al.l . 

Il998l : Iflebraud and Lootensl . [2005: Maranza no and Wagnerl. [ioOlal : [Fan" et all [20081 ). Visible shear in and above the 
shear thickening regime shows that shear thickening is not jamming in the sense of being associated with a yield 
stress or static structures (Figs. [6l [71 112|) . Below the onset of shear thickening, the particles may be settled or stuck 
together by attractions; in either case this corresponds to a locally static structure, and in many cases the system is 
jammed with a yield stress below the shear thickening regime. The shear thickening regime can be where more of 
the particles become involved in shear, i.e. becoming unjammed in the above sense, and it is the resulting dilation 
against the boundaries under shear which is responsible for Discontinuous Shear Thickening. 

One connection between Discontinuous Shear T hickening and jamming co mes from the observation that they are 
controlled by the same critical packing fraction (f> c (|Brown and Jaegerl . [2009). The shear rate at the onset of shear 
thickening goes to zero in the limit of <j> c , suggesting the limiting case of sh ear thickening correspond s to a yield stress, 
i.e. a jammed state, which is also what is found on the other side of C (I Brown and Ja eger. 200 j). It is suggestiv e 
that this is also near the critical point for the onset of dilation in dry granular systems (jOnoda and Linigerl . 11990). 
and the increase in dilation with shear rate must become steeper as this critical point is approached for our model to 
fit the data. Whether this feature is relevant for understanding dilation at non-zero shear rates in dense suspensions is 
uncle ar at this point . The other conn ection between Dis c ontinu ous Shear Thickening and jamming, which was hinted 
at bv lHolmes et alj (|2003l . [20051 ) and iMelrose and Ball |2004bl ). comes from the stress response, whose scales r max 
and Tj , respectively, are set by a confining stress due to the penetration of particles through the liquid-air interface 
in each case and which is transmitted through the system via force chains. This occ urs with dilation un der shear for 
Discontinuous Shear Thickening (Figs. [121 113p . and at rest for jammed suspensions (I Brown et al.l . [201ll ). In terms of 
the stress response, Discontinuous Shear Thickening could be considered to be a dynamic extension of jamming, but 
so far there is not yet a formalism for describing this. 

11.4 Summary of the mechanism for Discontinuous Shear Thickening 

The requirements for Discontinuous Shear Thickening under the mechanism we proposed, not necessarily limited to 
suspensions, can be summarized as: 

1. Dilation : The particle packing must attempt to dilate with increasing shear rate. 

2. Frustration : There must be geometric constraints that at least partially frustrate dilation. These constraints 
can come from surface tension at a liquid-air interface, the walls surrounding the system, or the particle stiffness, 
for example. This frustration usually requires a high packing fraction, close to the jamming transition. 

3. Confining stress: These constraints must provide a confining stress such that the most compliant constraint 
produces a restoring force that increases with dilational strain. 

4. Dominance : This confining stress must significantly exceed all stresses that prevent shear between grains and 
dilation, such as interparticle interactions or gravity. Otherwise, there is not enough stress increase from 
dilation to result in a positive slope on a viscosity curve and the global rheology may be shear thinning instead. 

Here we compare this mechanism to others that have been proposed for colloids and suspensions. Hof fman] (| 19821 ) 
argued that "[Discontinuous Shear Thickening] will occur in concentrated suspensions whenever the particles can 
segregate into layers parallel to planes of constant shear, but ar e constrained fro m free rotation at levels below some 
critical level of stress." Our model is in agreement with that of iHoffmanl (| 19821 ) in the sense that the onse t stress is 
determ ined by a point where the shear stress is large enough to shear particles in such a way to cause dilation. [Hoffman 
(|1982| ) argued for the same onset scaling law with the stress scale corresponding to the interaction between zeta 
po tentials of neigh boring particles. Our proposed scaling for the confining stress T max also gives the scale measured 
by IHoffmanl (|l972l ). whose data is included in Fig. [T5] While Hoffman observed an order- disorder transition along 
with Discontinuous Shear Thickening, it was shown that such a transition is not necessary jMaranzano and Wagnerl 
120021 : lEgres and Wagnerl . 120051 : lEgres et all 120061 ). which suggests the order-disorder transition is not the key to 

explaining the stresses. 

M any other papers have supp orted a hydrocluster mechanism (|Bradv and Bossisl . 1 19851 : iMaranzano and Wagnerl . 

l2001al : IShenov and Wagner. 2005) based on the ability of those models to predict some of the scalings for the onset 
Tmin- While those models have been successful at predicting weaker, continuous shear thickening, we suggest that 
for the much higher stress scales found in Discontinuous Shear Thickening to be reached, an additional mechanism is 
likely required to introduce another stress scale; for example, compression against the boundary introduces a confining 
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stress scale from the boundary stiffness. The hydrocluster model has been successful at describing the onset stress 
Tmin for Discontinuous Shear Thickening in some cases, but we have argued that the scalings for both r m ; n and 
Tmax predicted by different models for shear thickening can be more simply understood in terms of dominant stress 
scales, which are not dependent on whether stress transfer is through frictional or hydrodynamic interactions, or 
what mechanism initiates shear thickening. In this discussion, we have suggested that the dilational model based on 
suspensions could also apply to the colloid regime. While we have been able to compare our model to known scalings 
for T m i„ and to measurements of T max (Fig- 115) 1. more data is needed to directly test the dilational model in the 
colloid regime. 

12 Conclusions 

In this paper we proposed and experimentally validated a model for Discontinuous Shear Thickening in suspensions, 
which we briefly summarize here. The onset of shear thickening T m in occurs when the grain packing starts to dilate 
fFigs. rmil2p , which requires that the applied shear stress overcome any interparticle stresses preventing compressive 
shear between particles (Figs. [5j [6] [TJl ■ Dilation in turn causes the particles to push against the boundary, typically 
the liquid-air interface for suspensions open to the air, which pushes back with a restoring force to produce a confining 
stress on the suspension (Figs. 1131 1141 115)1 . The resulting normal stresses are transmitted through the packing via 
frictional interactions (Figs. [8] I9I10[) . resulting in a rapid increase in shear stress with shear rate corresponding to 
Discontinuous Shear Thickening. We generalized this shear thickening mechanism to other sources of a confining 
stress by showing that, when instead grains are confined by solid walls and have no liquid-air interface, T max is set 
by the stiffness of the most compliant boundary (Figs. [161 1171 [18)| . 

These results allowed us to suggest a generalization of the scaling laws for the stresses that bound the shear 
thickening regime to systems in other parameter regimes: the stress required for compressional shear between particles 
for T m i n (Sec. 111.2)1 . and the confining stress for T ma x (Sec. 111. 3) . We find that the experimental results can be 
described by a non-local constitutive relation where the stress does not come directly from the local shear rate, but 
where Discontinuous Shear Thickening comes out of a frictional term from the confining stress at the boundary which 
depends on the global dilation (Sec. Ill.Tj) . 
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14 Appendix: calculation of confining stress at liquid-air interface 
from dilation 

To connect the confining stress at the liquid- air interface to the measured dilation 5, a model is needed for the 
interface geometry. For simplicity we will assume spherical particles with a contact angle 6 = at the liquid-solid-air 
contact line since to obtain shear thickening the liquid must wet the particles. For this contact angle the interface 
geometry is equivalent to a sphere of air of radius r in contact with particles at the surface. We will calculate the 
geometry for a characteristic radius of curvature as if it is the same at each interstice between particles, and use 
these mean single-particle calculations as an estimate for the surface as a whole, ignoring variations in the surface 
curvature. The point of contact is defined by an angle a relative to horizontal as shown in Fig. 1191 The packing 
fraction of particles on the two dimensional surface will be represented by <j>2D- As an estimate we will use the value 
4>2D = 0.84 which corresponds to random close packing in two dimensions l|0'Hern et al.l . [2002h . This geometry gives 
enough constraints to relate the radius of curvature r to the contact angle. The dilation <5 can be connected to this 
geometry using conservation of volume. The confining stress scale from surface tension can then be calculated as 
proportional to 7/r as a function of S. 

From the geometry in Fig. 1191 the vertical components of the dimensions can be used to relate the radius of 
curvature to the contact angle 

r sin a — — (1 — sin a) + r m ■ (9) 
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Figure 19: Dimensions used for calculation of relationship between dilation 6 and confining stress due to 
surface tension at a liquid- air interface with radius of curvature r. 



The minimum radius of curvature r m comes from the minimu m rad ius of the insterstitial gap between particles. At 
an interstice between 3 particles in contact, this is r m /a = Wl/3 — 1/2 ~ 0.26 in three dimensions, while in two 
dimensions r m would be zero. The mean confining stress is modeled as 



7 / 2 
Tconf ~ — (1 — 02D Sin q) 



(10) 



where 1 — 4>2d sin 2 a is the fractional cross-sectional area around the outer edge covered by the liquid-air interface. 

Next, we relate the dilation 5 to a. Conservation of volume requires that the dilation match the enclosed volume 
of air per particle AV up to the maximum penetration of the particle when using the initial condition that the surface 
is flat (r = oo ) when 5 = 0: 



4</> 2D AV(r, a) 



The volume AV can be calculated as 



(11) 



AV = 2V cap + Vi (12) 
where V cap is the volume of a spherical cap interior to the point of contact (of which there are 2 per particle) 

„3 



V c , 



— cos Q ) 2 (2 + cos a) 



(13) 



and Vi comes from integrating the fluid volume in the mean surface normal direction from the furthest point of 
penetration of the particle up to the point of contact: 



1 

t>2D 



da 



1 — cos a cos(3a) — 9 cos a + 8 



4>2D 
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(14) 



With Eqns. [9l 1 1 01 and [TT] relating r con /, S, r, and a, they can be solved numerically to obtain a relationship between 
S and r con f, with the result shown in Fig. 1131 

These equations may be valid until the point where the contact line reaches the 2nd layer of particles from the 
surface, beyond which extra contact lines are made. This starts when the maximum penetration, defined in Fig. [19] 
as D — (a/2 + r) (1 — cos a) , reaches V3a/2 which corresponds to the layer width for a hexagonal packing. This limit 
is shown as the dashed line in Fig. 1131 
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